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ABSTRACT 


The  res-ssrch  described  has  been  directed  toward  an  understanding 
of  the  mechanics  of  rain  erosion  of  materials  traveling  at  very  high 
velocities.  An  analysis  of  the  motion  of  raindrops  in  a  hypersonic 
shock  layer  is  described.  A  high-speed  jet  impingement  analogy  of  The 
impact  of  liquid  drops  on  solid  surfaces  is  used  to  analytically  deter¬ 
mine  the  cavity  profile.  The  effect  of  stress  waves  generated  by 
liquid  impact  is  discussed.  It  is  shown  that  interned  damage  and 
spallation  resulting  from  stress  waves  is  often  much  greater  than  the 
material  failure  in  the  impact  area  itself.  Experiments  have  demon¬ 
strated  that  water-drop  impact  can  be  closely  simulated  by  the  use  of 
a  high-speed  water  jet.  Craters  formed  in  the  laboratory  by  the  jet 
are  compared  with  those  of  a  radome  surface  after  a  sled  test  in  a 


rainfield. 
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se:thx.  i. 

INTRt  AUCTION 

The  requirement  of  an  all-weather  capability  for  high-speed  air¬ 
craft  and  missiles  has  made  rain  ez-osior  protection  of  these  systems 
a  vital  design  factor.  No  dependable  design  criterion,  however,  has 
yet  been  formulated.  In  a  state-of-the-art  survey  of  raindrop  erosion, 
Wheelahan  (1967)  finds  that  very  little  quantitative  information  is 
available  on  damage  caused  by  raindrops  impacting  at  supersonic  veloc¬ 
ities.  Heymann  (1967)  says  that  there  exists  no  theory  for  predicting 
damage  to  be  expected  under  giver,  service  conditions.  A  -urge  number 
of  papers  and  reports  on  this  subject  has  been  published  u.cr  the  past 
few  years  and  many  of  them  have  contributed  toward  a  better  qualitative 
understanding  of  the  mechanism  of  rain  erosion  and  its  relationship  to 
material  properties,  although  most  of  them  have  been  limited  to  veloc¬ 
ities  less  than  1000  ft/sec.  These  papers  have  included  a  great  deal 
of  experimental  data  and  considerable  theoretical  work.  Many  of  the 
theories,  however,  have  been  rigorously  deduced  from  rather  arbitrary 
assumptions.  Contradictory  empirical  relations  have  been  formulated. 
One  only  needs  to  read  the  discussions  following  many  of  the  papers 
published  in  rhe  proceedings  of  various  rain  erosion  conferences  to 
see  that  there  is  much  disagreement  among  those  who  are  working  on  this 
problem.  There  is  considerable  disagreement  as  to  the  forces  that  are 
actually  responsible  for  the  erosion  and  cratering  resulting  from  a 
.liquid-solid  interaction.  According  to  Eisenberg  (1969),  "There  re¬ 
mains  the  mystery  of  which  mechanism,  shock  waves  or  jet  formation, 
is  responsible  for  actual  damage."  Engle  (1957),  Fyall  (1970),  and 
others  believe  that  the  damage  is  produced  by  the  introduction  of 
shear  and  tensile  stresses  as  a  consequence  of  the  high  locaiized  com¬ 
pressive  forces  exerted  by  the  liquid. 
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Regimes  of  Impact  Velocity 

The  various  theories  may  not  be  as  contradictory  as  they  first 
appear.  In  the  case  of  solid-solid  impact,  there  are  at  least  three, 
and  perhaps  four,  regimes  of  the  impact  phenomena  as  the  projectile 
velocity  increases.  These  are  illustrated  by  the  rather  idealized 

curve  Figure  1.1  which  shows  the  effect  of  velocity  upon  the  penetrat- 

* 

ing  power  of  projectiles  into  a  ductile  target.  In  the  1 '*w~velocity 
regime  the  penetration  is  strongly  influenced  by  the  strength  of  the 
two  materials.  In  this  regime  the  penetration  appears  to  increase  as 
the  4/3  power  of  the  velocity.  A  negligible  amount  of  the  projectile 
energy  is  converted  into  stress  waves.  A«  the  projectile  velocity  is 
increased  there  is  greater  deformation  and  fragmentation  of  the  pro¬ 
jectile.  At  some  point  there  is  no  increase  in  crater  depth  with  ir - 
creased  velocity;  in  fact  the  depth  actually  decreases  although  the 
crater  volume  may  continue  to  increase.  There  comes  a  point  wn^re  the 
depth  begins  to  again  increase  with  velocity.  The  crater  becomes  more 
hemispherical  in  shape  and  the  penetration  is  now  roughly  proportional 
to  the  2/3  power  of  the  velocity.  It  is  clear  that  a  transition  regime 
has  been  crossed.  The  new  regime  is  usually  called  the  hypervelocity- 
impact  stage.  '.ie  most  common  criterion  for  the  beginning  of  this  re¬ 
gime  is  when  the  ratio  of  projectile  velocity  to  that  of  a  stress  wave 
in  the  target  material  exceeds  unity.  At  these  velocities  the  stresses 
developed  by  the  projectile’s  deceleration  is  much  greater  tr.^r,  the 
yield  strength  of  either  the  projectile  or  target  and  the  materials' 
densities  affect  the  cratering  process  more  than  that  of  its  other 
properties  (mechanical  properties  such  as  target  and  projectile  strengths). 
Tr  the  velocity  is  further  increased,  an  impact  regime  is  reached  that  is. 
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as  yet,  relatively  unexplored.  It  is  called  the  regime  of  impact  explo¬ 
sions,  for  the  tremendous  amount  of  energy  becomes  sufficient  to  vaporize 
both  the  projectile  and  a  small  volume  of  the  target.  The  partitioning 
of  the  f.j.ojectile's  kinetic  energy  among  such  items  as  mechanical  defor¬ 
mation,  heat  generated,  stress  waves,  ejected  material,  and  perhaps  vapor¬ 
ization;  the  influence  of  material  properties ;  the  effects  of  stress  waves ; 
and  the  mechanism  of  penetration  and  cratering  are  entirely  different  for 
+he  various  impact  regimes.  It  is  apparent  that  there  are  also  different 
mechanisms  of  erosion  and  cratering  for  liquid-solid  impacts  at  various 
velocities.  Wheelahan  (19b7)  points  that  in  the  velocity  range  of 
1000  to  3000  ft/sec,  erosion  rather  than  cratering  is  the  dominant  form 
of  damage,  but  as  the  impact  velocity  is  increased  above  a  few  Mach  num- 
bars  the  damage  is  primarily  by  cratering.  Somewhere  between  these  two 
types  of  damage  there  is  obviously  a  regime  of  transition.  Those  appar¬ 
ently  contradicting  models  mentioned  nay  not  be  contradictory  at  all  but 
each  may  apply  at  various  impact  velocities  and  for  different  materials. 

For  a  particular  target  material  the  velocity-damage  curves  would  be  very 
different  for  liquid  and  solid  projectiles  but  v.ould  probably  converge 
in  the  hypervelocity  range.  If  this  be  true,  it  will  be  possible  to  simu¬ 
late  the  higher  velocity  liquid  impacts  with  solid  projectiles.  This, 
would  px-wbably  be  at  a  higher  velocity  than  is  presently  of  interest  in 
the  study  of  rain  damage  to  high-speed  aircraft  and  missiles.  For  example, 
the  hypervelocity  regime  of  aluminum  starts  at  a  velocity  greater  than 
20,000  ft/sec. 

A.1  tncugh  much  research  has  gone  into  the  accumulation  of  information 
now  available  relative  to  high-speed  liquid  impact,  it  is  obvious  that 
there  is  much  more  to  be  done.  Mathematical  models  must  be  derived  that 


take  into  account  the  effects  of  the  various  material  properties  including 
hardness, density,  elastic  modulus,  yield  strengths,  fracture  strengths, 
and  Poisson's  ratio,  as  well  as  surface  conditions  and  porosity;  the  frac¬ 
ture  criteria  for  various  materials,  both  ductile  and  brittle;  the  effects 
of  drop  size  and  angle  of  impact;  the  effects  of  larget  temperature ;  the 
impact  pressures  and  radial  flow  distribution  and  duration;  the  effects 
of  cavitation;  and  the  momentum  and  energy  imparted  to  the  target.  When 
iL  is  realized  that  this  must  be  done  for  the  entire  range  of  velocities, 
including  the  various  regimes  previously  discussed,  it  is  apparent  that 
much  is  yet  to  be  learned  relative  to  liquid  impacts.  This  will  require 
the  exchange  of  ideas  among  all  investigators  in  this  field  of  research. 
There  will  continue  to  be  conflicting  theories  but  controversy  always 
stimulates  research. 

Scope  of  the  Investigation 

Although  the  rain  field  may  be  defined  in  terms  of  drop  size,  dis¬ 
tribution,  velocity,  direction,  and  duration,  it  cannot  be  assumed  that 
this  describes  the  conditions  between  the  shock  wave  and  the  target  sur¬ 
face,  i  hich  in  this  study  is  a  conical  nosecone.  It  has  been  suggested 
that  water  drops  passing  through  the  shock  wave  may  be  fragmented, 
slowed  down,  and  turned  before  impacting  the  target.  The  shattering 
or  fragmenting  will  not  have  an  effect  if  the  time  to  traverse  the  shock 
layer  is  less  than  the  breakup  time.  If  the  deflectior  angles  are  of 
significant  magnitude,  the  total  number  of  impacts  will  be  reduced. 

If  the  raindrops  have  sufficient  momentum  to  prevent  them  from  being 
deflected,  all  raindrops  in  the  path  of  the  nosecone  will  impact.  Sec¬ 
tion  II  of  this  report  5s  an  analysis  of  the  motion  of  raindrops  in 
conical  shock  layers. 
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Two  damage-producing  attributes  of  liquid-drop  impact  are:  (a)  the 
impact  pressure,  and  (b)  the  sheet  jet  produced  by  the  radial  liquid  flow. 
Section  III  deals  with  the  mechanisms  of  cratering  and  erosion  resulting 
from  these  forces. 

A  third  source  of  damage  resulting  from  impact  is  the  stress  wave. 
This  wave  may  contribute  to  fractures  at  or  near  the  point  of  impact,  or 
may  produce  damage  by  reflection  from  the  target's  rear  surface.  These 
fractures  produced  by  reflected  stress  waves  will  generally  be  some  dis¬ 
tance  from  the  point  of  impact  and  when  they  occur  they  may  be  more  dam¬ 
aging  than  material  failure  in  or  near  the  impact  area.  This  is  an  im¬ 
portant  problem  as  it  may  dictate  the  required  nosecone  thickness.  It 
is  analyzed  in  Section  IV. 

This  investigation  includes  only  a  small  amount  of  experimental 
work.  A  water  jet  accelerator  was  constructed  and  targets  of  Lucite 
and  slip-cast  fused  silica  were  impacted.  A  description  of  the  equipment 
and  an  analysis  of  the  craters  and  fractures  constitute  Section  V. 

Finally,  an  SCFS  nosecone  was  examined  after  travel  through  a  rain- 
field  at  Holloman  Air  Force  Base  and  a  description  of  the  resulting 
damage  is  included  in  Section  VI. 

A  brief  description  of  proposed  future  work  is  given  in  Section  VII. 

As  the  various  sections  of  this  report  were  prepared  by  different 
individuals,  some  variations  in  symbols  and  forms  of  equations  will  be 
found;  hdwevef,  the  meanings  of  all  symbols  are  defined  in  the  text  and 
no  confusion  on  the  part  of  the  reader  should  result. 
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MOTION  OF  RAINDROPS  IN  A  HYPERSONIC  SHOCK  LAYER 

Introduction 

In  order  to  estimate  the  rain-erosion  damage  that  will  occur  when 
a  vehicle  flies  through  a  rain  field  it  is  first  necessary  to  knew  the 
relative  velocities  with  which  the  rain  drops  strike  the  vehicle  surface 
and  the  distribution  of  drops  along  the  surface.  The  theoretical  determin¬ 
ation  of  these  quantities  is  the  subject  of  this  section.  In  order  to 
compute  these  quantities  it  is  assumed  in  the  present  work  that  the  cloud 
of  drops  can  be  treated  as  a  quasi-continuous  gas-like  material.  This 
model  is  accurate  when  a  representative  volume  element  of  the  air-drop 
mixture,  small  compared  to  the  characteristic  dimensions  of  the  flow 
field,  contains  a  sufficient  number  of  drops  to  allow  meaningful  averages 
of  the  drop  properties  to  be  formed  over  the  volume  element.  These  aver¬ 
ages  are  then  treated  macroscopically  as  continuous  variables  and  the 
small  volume  is  assumed  to  have  the  mathematical  properties  of  a  differ¬ 
ential  element.  Thus  the  air-drop  mixture  is  viewed,  for  the  purposes  of 
macroscopic  analysis,  as  a  mixture  of  two  continuous  fluid-like  materials. 
General  discussions  of  the  governing  equations  appropriate  to  such  two- 
phase  flows  are  given  in  the  books  by  Wallis  (1969)  and  Soo  (1967)  and  the 
review  article  by  Marble  (1970). 

In  the  present  work  several  assumptions  have  been  used  to  develop  the 
governing  equations  to  be  discussed  subsequently.  These  are  the  subject 
of  this  paragraph.  First,  it  is  assumed  that  the  presence  of  the  drops 
has  an  insignificant  effect  on  the  motion  of  the  gas  phase.  This  assump¬ 
tion  is  justified  when  the  mass  fraction  of  the  dispersed  phase  is  small. 

Thu-  section  prepared  by  Dr.  John  Peddieson,  Jr.,  Assistant  Professor  of 
Engineering  Science  and  Chung-Hsein  Lyu,  Graduate  Instructor. 


This  condition  is  invariably  satisfied  in  shock  layers.  Second,  it  is 
assumed  that  the  volume  fraction  of  the  drop  phase  is  small  enough  so 
that  the  interphase  mass,  momentum,  and  heat  transfer  laws  can  be  eval¬ 
uated  using  formulas  appropriate  for  a  single  particle.  This  require¬ 
ment  is  more  stringent  than  the  first  one  but  it  appears  to  be  satis¬ 
fied  in  most  shock  layers.  Third,  it  is  assumed  that  collisicns  between 
drops  are  sufficiently  rare  so  that  the  drop-phase  stress  tensor  can  be 
ignored.  Finally,  it  is  assumed  that  the  drops  can  be  modeled,  with 
sufficient  accuracy,  as  identical  spheres  characterized  by  a  single 
effective  mass  and  radius,  both  of  which  are  independent  of  time.  Thus 
mass  loss  by  the  drops  is  not  explicitly  accounted  for.  Hass  loss  can 
occur  by  two  mechanisms:  evaporation  and  breakup.  In  a  state-of-the- 
art  survey  Wheelahan  (1967)  stated  that  neither  of  the  mechanisms  ap¬ 
peared  to  be  important.  A  set  of  equations  including  evaporation  ef¬ 
fects  is  given  by  Marble  (1970)  (  see  also  Panton  (1968),  Panton  and 
Oppenheim  (1968),  and  Lu  and  Chiu  (1966)).  An  order-of-magnitude  anal¬ 
ysis  of  these  equations  was  carried  out  in  the  course  of  the  Dresent 
work  and  tended  to  confirm  the  belief  that  loss  of  mass  due  to  evapora¬ 
tion  is  negligible.  Recent  experimental  research  on  drop  breakup  (see, 
for  instance,  Waldeman,  Reinecke,  and  Glenn  (1972))  indicates  that  the 
mass  loss  associated  with  this  phenomenon  can  be  appreciable  at  high 
Mach  numbers.  Unfortunately  the  existing  experimental  data  has  not  yet 
been  put  in  a  form  suitable  for  incorporation  in  the  theory  to  be  subse¬ 
quently  discussed.  An  expression  relating  the  rate  at  which  mass  is  lost 
by  a  single  drop  to  the  parameters  of  the  flow  of  the  gas  relative  to  the 
drop  is  required.  To  date,  no  such  relationship  has  appeared  in  the 
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Uterature.  Thus  drop  breakup  cannot  be  considered  explicitly.  A 
possible  way  of  correcting  the  theory  for  drop  breakup  will  be  dis¬ 
cussed  later. 

Previous  investigations  along  the  lines  of  the  present  work  are 
United  to  the  following  papers.  Probstein  and  Fassio  (1970)  computed 
the  drop-phase  velocity  Held  and  streanline  pattern  for  shock  layers 
on  thin  wedges  and  cones  and  in  the  stagnation  regions  of  cylinders 
and  spheres.  They  employed  the  constant  density  appraxination  to  deter¬ 
mine  the  gas  flow  field.  For  the  wedge  they  obtained  results  in  closed 
font  while  numerical  integration  was  necessary  to  solve  the  equations  in 
the  other  cases,  tfaldeman  and  Reinecke  (1971)  obtained  results  for  the 
velocity  field  and  streamline  pattern  for  flow  past  a  cone  and  a  sphere. 
They  relaxed  the  constant  density  assumption  invoked  by  Probstein  and 
Fassio  (1970)  but  solved  the  governing  equations  for  the  drop  phase 
approximately  by  a  perturbation  method.  Both  of  these  reports  also 
gave  results  for  the  collection  efficiencies  of  the  various  body  shapes 
analyzed  in  them.  (The  collection  efficiency  is  defined  to  be  the  ratio 
of  the  actual  rate  at  which  aasr  is  removed  from  the  drop  phase  by  col¬ 
lection  of  drops  on  the  body  surface  to  the  rate  that  would  exist  if  the 
motion  of  the  drop  phase  were  entirely  unaffected  by  the  presence  of  the 
shock  layer. )  Further  collection  efficiency  calculations  were  made  by 
Spurk  and  Gerber  (1972)  who  considered  cones  and  slender  power- law  axi- 
symmetric  bodies. 

In  the  present  work  the  drop-phase  velocity  field,  strec-,1 ine  pat¬ 
tern,  density  distribution,  and  temperature  distribution  are  '.omputed 
for  conical  shock  layers.  Since  the  collection  efficiency  seems  to  have 
little  relevance  to  radome  design,  no  results  for  this  quantity  are 
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presented.  These  calculations  should  be  representative  of  tfct  condi¬ 
tions  occurring  in  radcae  shock  layers  because  the  results  of  Spurk 
and  Gerber  (1972)  show  that  snail  deviations  fron  a  perfectly  conical 
body  shape  produce  insignificant  changes  in  the  floe  variables  of  the 
dispersed  phase.  The  solutions  are  obtained  nunericaily  and  results  are 
presented  graphically  to  illustrate  various  paranetric  trends.  In  addi¬ 
tion,  sone  results  are  given  for  shock  layers  an  wedges.  While  such 
calculations  have  no  direct  application  to  radone  technology,  the  wedge 
is  a  useful  test  case  because  the  gas-phase  flow  field  has  an  especially 
sinple  fore.  In  fact,  by  invoking  the  thin-shock- layer  approximation  it 
is  possible  to  obtain  solutions  in  closed  fbm  for  this  case.  This  pro¬ 
cess  was  began  by  Probstein  and  Fassio  (1970)  and  is  completed  in  the 
present  work. 


Governing  Equations 


The  appropriate  governing  equations  for  analysis  of  the  behavior 
of  the  drop  phase  are  found  from  balances  of  mass,  linear  momentum,  and 
energy.  These  are  respectively 
pP>t  *  = 0 

pp{Vt +  v^V =  ? 
pp(Vt +  W =  s 

where  (2.1a)  results  from  balance  of  mass,  (2.1b)  from  balance  of  linear 


(2.1) 


momentum,  and  (2.1c)  fron  balance  of  energy;  is  the  drop-phase  mass 
density,  is  the  drop-phase  velocity  vector,  ep  i3  the  drop-phase 
internal  energy,  f  is  the  force  per  unit  volume  applied  to  the  drop  phase 


by  the  gas  phase,  S  is  the  heat  per  unit  volume  transferred  to  the  drop 
phase  by  the  gas  phase,  V  is  the  gradient  operator,  t  is  time,  and  a 


coma  denotes  partial  differentiation  villi  respect  to  the  Miwfag 
«4»c*apt.  To  reader  (2.1)  determinate  oanstittitiv*  equations  for  y 
P»  «d  S  »«  required.  For  Hie  internal  enemy  it  is  aesaeed  (see 
Sarfcle  (197?))  that 

'p  *  cPTp  OM 

wfw  Cp  is  the  heat  capacity  of  the  dispersed  phase  and  Tp  is  its  tea- 

jerature.  Tito  sets  of  acatntus  and  heat  supply  terns  are  used  in  Hie 

present  work.  The  first  is  (vith  u  denoting  Hie  fluid-phase  velocity  vector) 

?  ^  -(PpA1)|up-aj1*1>(Up-n) 

S  =  (pc/i,)(T  -T)  (2.3) 

P  2  p 

wb e.-e  c  is  the  specific  heat  at  constant  pressure  of  the  gas  please, 
b  is  a  numerical  factoi;  and  and  t2  are  aonentus  and  tberaal  relaxation 
factor.'  defined  as 
ij  =  2^*‘ x/Jlia^p^V 

t2  =  nc/4»a*  (2.4) 

In  (2.4)  N  is  »  ruverical  factor,  a  is  the  mass  and  a  is  the  vadim  of 
a  representative  spLerical  particle,  p  is  the  gas-phase  density,  p  is 
the  gas-phase  viscosity,  and  e  It  the  gas-phase  theraal  conductivity. 

The  relative  Reynolds  mn'ier  is 

R  =  (2ap/u)jup-uj  (2.5) 

Its  value  determines  the  numerical  factors  b  and  H  as  follows. 

0  <  R  <  1  :  b  *  1,  M  *  24 

1  <  R  <  1000  :  b  =  0.6,  N  =  24 

1000  <  R  <  **  :  b  =  0,  N  =  0.44  (2.6) 

Equation  (2.3a)  is  obtained  by  dividing  the  standard  drag  curve  for 

steady  flow  past  a  sphere  into  three  parti  end  fitting  an  equation  to 
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eadi  part  as  described  by  Pnefcsteia  and  Fassio  (i«70).  It  JU  erased 
in  the  present  «oric  that  tbe  only  eifdfiat  contribution  to  tb*  inter- 
phase  force  is  tbe  steady  state  drag.  This  is  knows  to  predeninate  the 
drag  doe  to  added  aass  and  that  doe  to  instead?  notion  in  aost  cases, 
lift  farces  were  included  in  sane  of  the  calculations  nade  in  tbe  course 
of  the  present  work  and  were  found  to  bane  a  negligible  effect-  Equa¬ 
tion  (2 -3b)  results  frow  the  pure-condoction  beat  transfer  law  for  a 
single  sphere.  Equations  (2.3)  have  a  convenient  analytic*’*  farm  which 
facilitated  the  deteminatica  of  the  exact  solution  referred  to  previously. 

These  equations  were  useful  in  checking  the  conputer  prograw  developed  in 
the  present  work. 

The  second  set  of  drag  and  heat  transfer  laws  used  is 
P  =  -(Pp/XjMl  ♦  (64/R^5 ) ) j Up-u|  (up-u) 

S  =  -(ppc/T2)(l  +(3/10)R1/2Pr1/3)(Tp-T)  (2.7) 

where  all  sysbols  previously  defined  retain  their  neanings  with  the  stipu¬ 
lation  that  b  be  set  equal  to  zero  in  (2.4a).  In  (2.7b)  Pr  is  the  Prandtl 
□uaber  defined  by 

Pr  =  uc/r  (2.8) 

foraulas  of  tbe  font  of  (2.7)  have  been  proposed  by  nany  authors  (see,  for 
instance.  Panton  (1968))  with  slightly  different  values  of  the  coefficients 
and  exponents  and  are  known  to  be  reasonably  accurate  over  a  wide  range  of 
Reynolds  and  Prandtl  numbers . 

Thin  Shock  Layers  on  Wedges* 

Consider  a  wedge  of  surface  length  h  and  half-angle  6fa.  The  free¬ 
st  ream  velocity  is  and  the  ratio  of  the  gas-phase  density  before  the 
shock  to  that  a^ter  the  shock  is  e.  The  density  of  the  gas  phase  and  the 


*A  shortened  account  of  this  work  is  contained  in  Peddieson  and  Lyu  (1973). 
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drop  phase  just  Mind  the  shock  are  pg  and  respectively.  The 
free-strsae  tespcrature  of  tike  Mixture  is  T—.  It  is  convaafent  to 


to  use 


cartesian  coordinates  x  and  y  with  x 


parallel  to  the  body 


face  and  y 


appeal  to  it.  It  is  convenient  to  define  dinension— 


variables  as  follows. 


s  =  x/L,  n  =  y/L 

Fes  +  ^  =  u/0-»  rp®»  *  Gpen  =  "p70- 
Q  *  e/P8,  Qp  *  Pp/Pp. 

h  =  t/t.,  Hp  =  t/t. 

where  t  and  ?  are  unit  vectors  associated  with  the  --  and  n  axes, 
s  n 

Substituting  (2.9)  and  (2.10)  into  (2.1),  (2.2),  ud  (2.3)  (equations 
(2.3)  are  used  initially  because  their  fauns  facilitate  the  subsequent 
analytical  work),  and  substituting  (2.2)  and  (2.3)  into  (2.1)  results  in 

'Vp'-S  +  %%K  -  0 

f-F  ,  .  GF  .  «,((F-F  )2  +  (G-GTi)2)<1_b,/2(F-F  ) 

p  p  s  p  pn  1  P  p  p 

F  G  ,  +  6  G  »„  =  «i<(r~Fn>2  +  (G-Gn)2)(1"b>/2(G-G  ) 

ppsppnl  p  p  p 

FH  .  ♦  GH,  =  o~(H-H  ) 

p  p’s  p  p’n  a  p 

where 

=  L/t^,  oj  =  (L/xoUj(c/ep) 

The  gas  flow  variables  are  modeled  by  the  inviscid  solution  for  super¬ 
sonic  flow  past  a  wedge.  This  solution  is  known  exactly  and  for  a  thin 
shock  layer  (which  oceurs  when  6^  and  c  are  snail)  ?..as  the  forn 

F  =  1,  G  =  0,  H=H,  Q  =  1 

s 

where  the  subscript  s  denotes  the  value  of  the  variable  on  the  downstream 
side  of  the  shock  wave.  The  shock  wave  is  known  to  be  a  straight  line 


(2.9) 


(2-10) 


(2.1J ) 


(2.12) 


(2.13) 


making  an  angle  (for  a  thin  shock  layer) 


-12- 


%  =  v(1  - c)  < 

with  title  centerline  of  the  wedge.  The  solution  of  (2.11)  wost  satisfy 
the  conditions 

rP  =  1-  V'V  V1’  V1  ‘ 

on  the  shock  wave  which  is  the  line 

n  =  (es  -  ( 

This  problem  is  nost  easily  solved  using  a  Van  Rises  transformation 
fail  liar  fna  boundary-layer  theory.  This  is  done  by  defining  a  stream 
function  *p  such  that 

*p’s  '  *V*n  '  %Fp  * 

and  converting  the  independent  variables  from  s  n  to  f>  and  z  where 

z  =  n.  The  total  differential  of  +  is 

P 

d*p  =  ♦p*sds  *  V®**  s  *  Vp*1  ( 

Solving  (2.18)  for  ds  yields 

ds  =  -(l/QpGpJd^  ♦  (rp/Gp)dz  ( 

How  regarding  s  as  a  dependent  variable  which  is  a  function  of  and  z 
it  is  clear  that  .  ■ 

=  •»  =  -WpGp  •  S-2  *  VS  ( 

l 

Converting  the  derivatives  in  (2.11b,c,d)  and  rewriting  (2.20b)  one  obtains 
Vp-z  *  “i<(1-rp>2  *  Gp2)<1M/2(i-rp) 

vv*  *  -v^-v2  *  y>‘M)/% 

W*  "  «2<VV 


(2.1%) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


(2.20) 


Gp*»z  =  Fp 

where  (2.13a,b,c)  have  been  used.  Substituting  (2. 15a,bfd)  and  (2.16)  into 
(2.18),  integrating  along  the  shock  wave,  and  assuming  ^(0,0)  =  0  yields 
i  =  0_z_/(9_-8,  ) 

p  S  S  S  D 


(2.21) 


(2.22) 


-1%- 


txa  tiiis  Is  dene  (2.20a)  reads 

*  -V«VebV-*,) 

The  Initial  conditions  are 

V1-  w  V1*  * 1  wv 

on  *  =  Since  equations  (2.21)  contain  no  derivatives  with  respect 
to  Zg  they  cm  be  integrated  like  ordinary  differential  equations  along 
the  stneanlines  (lines  of  constant  x%).  The  dependence  of  the  solution 
an  Zg  enters  thrcugh  the  application  of  the  boundary  conditions.  For  all 
values  of  b  the  solution  far  the  tangential  velocity  is 

fp  =  1 

The  other  variables  have  different  fores  depending  an  the  value  of  b. 

They  are 

b  =  1  : 

s  =  (V<es-V>  "  (i/»i)i°g(i  -  («l/«b)(ss-*)) 

Gp  =  -®b  ♦  °l(2s-z) 

Hp  -  «s  -  (Hs-l)(l  -  (a1/eb)(zs-*))°2/al 

Qp  =  1/(1  -  <V®s>(1  "  («i/eb)(zs”2))) 
b  =  0.6  : 

s  =  (zs/(es-0b)>  -  (5/2a1©b2/^)(l  -  (1  -  (3^/5^  3/5)  (2s-z))"2/3) 

Gp  =  -^(1  -  (3o1/50b3/5)(zs-2))5/3 

Hp  =  «s  -  (Hs-l)siip(-(5i«2/2«1^)2/5)(a  -  (3o1/5eb3/5)(2s-z))"2/3-X)) 

Qp  =  1/(1  -  (ejj/esMl  -  (1  -  (3o1/56b3/5)(zs-z))5/3)) 
b  -  0  : 

s  =  <zs/<V0s))  ~  -  exp(a1(zs-z))) 

Gp  =  -0j,  expC-a^Zg-z)) 

Hp  =  Hg  -  (Hg  -  Dexpt-toj/cj^BjjXexpt^tZg-z))-!)) 

Qp  =  1/(1  -  (e^/egHl  -  exj>(d}(zs-z)))) 


(2.23) 

(2.2%) 

(2.25) 


(2.26) 

(2.27) 


(2.28) 


Setnlts  equivalent  to  (2.26ajb),  (2.27*,b),  and  (2.2tajb)  have  been  given 


by  Probstein  and  Fassio  (1970).  They  used  a  Lagrangian  description.  The 
temperature  and  density  results  are  new.  The  use  of  m  Eulerian  descrip¬ 
tion  facilitates  the  deterndnation  of  the  drop-phase  density.  These  re¬ 
sults  were  used  as  a  check  on  the  niawrical  calculations  to  be  described 
subsequently  and  a  check  on  the  apprqadnatc  analytical  solutions  discussed 
in  Appendix  1. 

Munerical  Solutions  for  Thin  Shock  Layers 

In  this  sub-section  the  synbol  j  is  introduced  to  differentiate  be¬ 
tween  plane  and  axisyetric  flow.  For  the  fomcr  j  =  0  while  for  the 
latter  j  =  1.  The  discussion  will  be  for  the  conical  sase  but  the  result¬ 
ing  equations  are  also  valid  for  wedge  flows  if  j  is  set  equal  to  zero. 

For  numerical  work  it  was  found  convenient  to  use  a  systew  of  polar 
coordinates  r  and  8  with  r  =  0  at  the  leading  edge  of  the  body  and  0  =  0 
on  the  axis  of  syunetry  of  the  cone.  Defining  the  diwens iooless  variables 
£  -  r/L,  n  =  8 

Fe  +  Ge  =  u/U  ,  F  e  +  G  e  =  u  /U 

where  e^  and  e^  are  unit  vectors  associated  with  £  and  n,  substituting 
(2.29)  and  (2.10)  into  (2.1),  (2.2),  and  (2.3)  and  substituting  (2.2)  and 
(2.3)  into  (2.1)  yields 
(f,<tn)J‘Qplp),(  .  ((Zn^QpGp)^  =  0 

VW  *  <Gp/5)Fp-r,  '  <Gp2/U  =  “l((F-Fp)2  *  (G-Gp)2)<1'b)/2<F-rp) 

W?  *  V{>Gp’P  ♦  <rpGp/5>  *  “l<(r-Fp>2  +  <G-G?)2)a-b)/2(G-Gp) 

W?  *  <Gp/5)«p,p  *  »2(HS-Hp) 

The  gas  flow  variables  are  modeled  by  the  constant  density  approximation 
(see  Probstein  and  Fassio  (1970)) 


(2.31) 


F  *  1,  G  =  -teOt-^))*,  H  -  Q*1 

(Bote  that,  within  the  thin-shock -layer  approximations,  these  egressions 
are  exact  for  j  *  0.)  The  shock  van  is  a  cone  making  an  angle 

•s  =  ej/d  -  (c/2^)) 

with  the  axis  of  syametxy.  The  conditions 

V1’  vv  V1*  %sl 

must  be  satisfied  on  the  shock  ware. 

Again  it  is  useful  to  employ  a  Von  Mises  transformation.  Defining 
a  stre.  fmcti™  *p  such  ttat 

V?  *  -«1>\V  Vn  = 

and  converting  to  the  independent  variables  t  and  z  =  n  -  6^  one  obtains 

GPFP*z  •  V  =  al€(<1  ‘  FP>2  +((2z)j+  Gp)2)(1_b)/2(i  *  Fp> 

GLG  ,  +  F  G  =  -a,£((l-F  )2  +((2zP+  G  )2)(l-b)/2((2z)5+G„) 
ppzpplp  p  p 

6PHP»Z  =  °2<Hs  -  V 

V’«  =  5FP 


(2.32) 


(2.33) 


(2.34) 


Op  =  -l/(U'z  + 


(2.35) 


The  relationship  between  f  and  £  the  shock  position  where  a  streamline 

P  s 

originates  is  easily  found  using  (2.34)  to  be 

*D  =  (0  e1+j/(i+j)K  1+j 

Using  this  expression  (2.35)  can  be  written  as 

Qp  =  -<{s/s)j<V<*  *  ’ 

s 

Equations  (2.35a,b,d)  were  numerically  integrated  simultaneously  along 
various  streamlines  using  a  fourth-order  Runge-Kutta  algorithm.  Since 


(2.36) 


(2.37) 


each  streamline  represents  a  line  of  £g  =  a  constant,  the  partial  deriva¬ 
tives  with  respect  to  z  can  be  treated  as  ordinary  derivatives  in  this 
calculation.  After  the  determination  of  Fp,  Gp,  and  £  (2.35c)  was  integrated 
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using  a  fourth-order  Runge-Kutta  method  to  obtain  Hp  and  Qp  was  obtained 
from  (2.37)  where  the  derivative  with  respect  to  Cg  was  approximated  by 
a  two-point  backward  difference.  A  discussion  of  the  numerical  results  • 
obtained  by  this  method  is  the  subject  of  the  next  three  paragraphs . 

Figures  2.1  and  2.2  show  the  distribution  of  the  drop-phase  normal 
velocity  and  the  drop-phase  density  along  the  vehicle  surface.  Ihese 
calculations  are  for  a  vertex  half-angle  0fa  =  0.1  and  a  gas  density  ratio 
e  =  0.1.  Results  are  plotted  for  various  values  of  the  dimensionless 
drag  parameter  a  .  These  results  are  for  the  linear  drag  law  b  =  1. 

Figure  2.1  pertains  to  wedges  while  Figure  2.2  pertains  to  couical  shock 
layers.  It  can  be  seen  that  for  small  values  of  the  drag  parameter  the 
drop-phase  normal  velocity  and  the  drop-phase  density  are  nearly  constant 
along  the  surface.  This  represents  a  situation  in  which  the  motion  of 
the  drops  is  almost  unaffected  by  the  presence  of  the  shock  layer.  As 
increases  significant  deviations  from  the  free-stream  values  of  the 
variables  quickly  become  apparent.  It  is  clear  that  this  effect  is 
greater  for  the  wedge  than  for  the  cone.  Figure  2.1  represents  either 
the  numerical  solution  or  the  analytical  solution  represented  by  (2.26a,d). 
For  the  step  sizes  used  in  this  work  (AC  =  .005,  Az  =  ((0  -0,  )/100))  it 

S  S  D 

was  found  that  the  numerical  and  analytical  solutions  were  indistinguish¬ 
able.  Reductions  in  the  step  sizes  failed  to  produce  any  change  in  the 
solution  variables. 

Figures  2.3  through  2.8  present  further  results  for  the  problem  dis¬ 
cussed  in  the  previous  paragraph.  The  distributions  of  various  variables 
are  plotted  along  selected  streamlines.  The  parameter  4>  =^/(9(l+j)) 
increases  as  the  distance  from  the  origin  of  the  streamline  to  the  nose 


increases.  Figures  2.3  and  2.4  show  the  variation  of  the  normal  velocity 
with  z  for  two  values  of  o^.  On  the  streamlines  close  to  the  nose  of  the 
vehicle  the  particles  are  slowed  down  the  least.  The  drag  becomes  in¬ 
creasingly  more  effective  as  the  distance  from  the  nose  increases.  This 
trend  becomes  more  pronounced  as  ot^  becomes  larger.  As  the  normal  veloc¬ 
ity  is  decreased  along  a  given  drop-phase  streamline  the  drop-phase  density 
increases.  This  is  shown  in  Figures  2.5  and  2.6.  The  increases  in  density 
are  larger  for  plane  flow  than  for  axisymmetric  flow  and  become  steadily 
larger  as  the  value  of  dimensionless  drag  coefficient  is  raised.  Figures 
2.7  and  2.8  depict  selected  drop-phase  streamlines.  The  streamlines  are 
slightly  deflected  downstream  by  the  action  of  the  interphase  drag.  This 
effect  is  accentuated  as  the  distance  of  the  origin  of  the  streamline  from 
the  nose  increases  and  as  a1  increases.  The  radial  velocity  Fp  has  not 
been  plotted  because  its  value  was  found  to  be  very  close  to  unity  for  all 
the  calculations  discussed  here. 

Figure  2.9  shows  some  comparisons  of  the  results  obtained  in  different 
drag  regimes  (different  values  of  b)  for  =  0.1  and  e  =  O.J.,  Obviously 
the  sensitivity  of  the  solution  variables  to  the  value  of  b  increases  as 
a-^  increases  because  b  is  contained  only  in  the  drag  terms  whose  magnitude 
is  determined  by  the  value  of  a^.  It  is  clear  from  the  plotted  data  that 
for  =  10  significantly  different  values  of  the  flow  variables  exist  in 
different  regions  of  the  standard  drag  curve.  For  thin  shock  layers  it  is 
reasonable  to  assume  that  one  drag  law  will  be  applicable  throughout  the 
flow  field.  This  has  been  verified  by  the  calculations  of  Spurk  and 
Gerber  (1972).  Calculations  of  the  type  described  in  the  previous  two 
paragraphs  were  carried  out  for  the  oases  of  b  =  0.6  and  b  =  0.  Since 
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the  results  were  qualitatively  siailar  to  those  already  discussed  they 
are  not  presented  graphically. 

Numerical  Solutions  for  General  Conical  Shock  Layers 

It  is  now  desired  to  eliminate  any  assuaption  that  either  the  cone 
or  the  shock  layer  is  thin.  Since  this  iaplies  that  a  wide  range  of 
relative  Reynolds  mnfcers  nay  be  encountered  by  a  typical  drop  as  it 
traverses  the  shock  layer,  it  is  appropriate  to  use  the  drag  and  heat 
transfer  laws  given  by  (2.10).  The  governing  equations  (2.30)  are 
now  replaced  by 

U2  sin(u)QpFp)tJf  ♦  (e  sin(n)QpFp),^  =  0 

FPFP’C  +  (GP/5)FP*n  '  <Gp/5)  =  °1W(1  +(64/(°3<3((f-fp)2  + 
(G-Gp)2)1/2)4/5))((F-Fp)2  +  (G-Gp)2)1/2(F-Fp) 

FPGP*S  +  <GP/5>Gp*n  +  <FpGp/0  S  °15Q<1  +<64/^3Q^F“Fp>2  + 

(G-Gp)2)1/2)4/5))((F-Fp)2  +  (G-Gp)2)1/2(G-Gp) 

FpHp,c  +  (Gp/UHp,n  =  o^U  +  (3Pr1/3/iq)P3Q((F-Fp)2  + 

(G-Gp)2)1/2)1/2)(H-Hp) 

where  is  given  by  (2.12)  with  b  =  0  and  p  replaced  by  pg  in  the  equa¬ 
tion  defining  (equation  (2.4a)  with  b  =  0)  and  a3  =  (2aUfl0ps)/p.  Note 
that  ag  has  the  form  of  a  Reynolds  number  based  on  the  free-stream  veloc¬ 
ity,  the  effective  particle  radius,  the  gas  density  just  behind  the  shock, 
and  the  gas  viscosity.  The  gas  flow  variables  are  modeled  as  follcws. 
f  =  cos(es)  -  ajQ  -  (e/es)) 

G  =  -a2(sin(03)/es)(e-eb) 

Q  =  1  +  a3(l  -  (0/0s)) 

H  *  Hs  +  a4(l  -  (0/6g ) ) 


(2.38) 


l 


(2.39) 
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Equations  (2.39a,b,c)  were  obtained  by  Waldaan  and  Reinecke  (1971)  by 
fitting  linear  functions  to  existing  theoretical  and  experimental  data 
for  the  gas  flow  field.  Equation  (2.39d)  was  obtained  the  sane  way  in 
the  course  of  the  present  work.  Representative  values  of  the  constants 
contained  in  (2.38)  are  a^  =  0,  a2  =  1.67,  ag  =  0.35,  84  =  4.0.  Calcula¬ 
tions  carried  out  in  the  course  of  this  work  showed  that  the  drop-phase 
flow  variables  were  rather  insensitive  to  the  values  of  the  a's.  The 
shock  angle  is  computed  (following  Waldman  and  Reinecke  (1971))  from  the 
formula 


es=eb/(l  -  (e/a2)) 

The  boundary  conditions  to  be  satisfied  on  the  shock  (0  =  0S)  are 

Fp  =  cos(%),  Gp  =  -sin (0S),  Hp  =  Hs,  Qp  =  1 

In  this  case  the  drop-phase  stream  function  defined  so  that 

Vs =  sin(n)QpV  Vu s  *>2  sin(n)QpFp 

and  the  modified  governing  equations ( 2 . 35 )  are  replaced  by 

G F  ,  -  G2  =  a.CQ(l  +  (64/(c.  Q((F-F)2  + 

P  P  z  p  1  3  P 

(G-G_)2)1/2)4/5))((F-F_)2  +  (G-G  )2)1/?(F-F  ) 

P  P  p  P 

Vy2  +  FPGp  =  VQ(1  +(64/(a3Q<(F-Fp)2  ♦ 

(G-Gp)2)1/2)4/5))((F-Fp)2  +  (G-Gp)2)1/2(G-Gp) 

GpHp,2  =  o2C(l  +  ( 31^/3/10)  (bt3Q((F-Fp)2  +  (G-Gp)2)1/2)1/2)(H-Hp) 

G  5.  =  CF 

p  z  p 

Qp  =  -1/(^ sin(z  +  W** } 

Equation  (2.36)  is  replaced  by 
*  =  ((sin(8  ))/2)£2 

P  00 

This  relationship  can  be  used  to  rewrite  (2.43e)  as 
Qp  =  -(Cs/5)(sin(es)/sin(z+eb))(sin(es)/GpS,£  ) 


(2.40) 


(2.41) 


(2.42) 


(2.43) 


(2.44) 


(2.45) 
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Equations  (2.43a,b,c,d)  and  equation  (2.45)  were  solved  numerically  using 

the  method  discussed  earlier.  These  solutions  were  found  to  be  insign  if  1- 

*  :  ! 
cantly  different  from  those  described  earlier,  for  corresponding  values  r* 

! 

the  parameters . _  i 

i  .  I 

In  order  to  carry  out  a  parametric  study  a  representative  density 
ratio  of  e  *  0.2  was  selected  and  used  in  all  the  calculations  to  be  eon- 

•  *  i 

■  • 
i 

sidered  below.  Figures  2.10  and  2.11  illustrate  the  influence  of  the 

I  !  • 

parameter  on  the  solution  variables  evaluated  on  the  wall  of  the 'cane 

for  6jj  =  0.1.  Far  large  a3  the  quadratic  drag  law  is  approached  and  the  : 

solution  is  insensitive  to  a3.  As  a3  decreases  the;  sensitivity  of  the 

■ 

solution  to  its  value  increases  and  large  reductions  in  the  drop-phase 
normal  velocity  and  increases  in  the  drop-phase  density  occur.  The  power 
of  to  influence  the  numerical  results  increases  with  increasing  c^. 

Figures  2.12  through  2.17  depict  the  variation  of  the  transverse 
velocity  and  the  density  along  selected  streamlines  and  the  shapes  of 

i 

these  streamlines  for  various  values  of  03.  The  vertex  half- angle  for* 
these  cases  is  0^  =  0.2.  The  parameter  =  £g/10.  It  is  evident  that 

i 

the  variation  of  quantities  from  streamline  to  streamline  is  increased 
by  decreasing  o3»  F<pr  thq  data  plotted  in  these  figures  there  is  little 
variation  across  the  streamlines  for  any  of  the  parametric  values  shown. 
Comparison  of  the  results  for  O3  =  105  and  a3  =  10  9  reveals  that  they, are 

virtually  identical  indicating  that  the  quadratic-drag- law  asymptote  to 

1 

the  standard  drag  curve  has  been  attained. 

Figures  2.18  through  2.20  show  the  influence  of  ^he  cone  half-angle 
Sjj  on  the  surface  values  of  the  normal  velocity  and  the  density  of  the 
drop  phase  for  various  values  of  03.  As  the  angle  0^  increases  the 
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amoupt  of  variation  of  the  flow  variables  aloof  the  surfaces  decreases. 
This  trend  is  aost  apparent  for  =  102.  It  is  obvious  that  both  the 
normal  velocity  and  the  density  of  the  drop  phase  are  increasing  functions 


of  tiie  cone  half-angle. 

Figures  2.21  through  2.23  illustrate  the  influence  of  t,  on  the  drop- 
phase  nornal  velocity  and  the  drop-phase  density  at  the  vehicle  surface 
for  6^  =  0.4  and  various  values  of  a^.  As  discussed  earlier  snail  values 
of  lead  to  virtually  constant  distributions  of  and  along  the 
surface.  As  increases  large  variations  in  the  values  of  the  flea* 
variables  at  the  wall  appear.  This  effect  is  nuch  none  pronounced  in 
this  case  than  it  was  for  the  thin  shock  layers  discussed  previously. 

The  largest  variations  occur  for  the  smallest  value  of  the  Reynolds 
number  parameter  <*3.  Even  for  the  larger  values  of  <*3  significant  varia¬ 
tions  of  and  Q  can  occur. 

P  p 

Figures  2.24  through  2.29  present  results  for  the  streamline  shapes 
and  the  variations  of  the  transverse  velocity  and  the  density  of  the  drop 
phase  along  selected  streamlines  for  the  same  geometry  discussed  in  the 
previous  paragraph.  The  streamline  parameter  i|<  =  £s/10.  In  all  cases  it 
can  be  seen  that  the  variations  of  the  plotted  variables  with  z  increases 
as  the  distance  of  the  streamline  from  the  leading  edge  of  the  cone  in¬ 
creases.  This  trend  becomes  less  prominent  as  is  increased  and  more 
prominent  as  is  increased. 

In  all  the  cases  discussed  in  the  present  work  the  temperature  dis¬ 
tributions  were  computed  along  with  the  other  results  that  have  been 
discussed.  Since  the  surface  temperature  distribution  of  the  drop  phase 
does  not  now  appear  to  be  pertinent  to  design  procedures  for  the  preven¬ 
tion  of  rain  erosion  damage  this  data  was  not  presented  graphically. 
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For  drop  sizes  typical  of  those  food  in  rain  fields  wpwatatiit 
nines  of  the  diwarfonless  drag  parsetcr  «j  =  L'r.  lie  is  the  range 
HT1  -  dj  -  1  if  tj  *  */(.22*p#a2)  *  6.0%aCcp/p#>  Cflp  is  the  density  of 
the  Material  contained  in  a  drop)  is  enlsiited  naiag  the  initial 

drop  radios.  It  should  be  possible  to  account  for  breakup  of  the  drops 
(when  this  occurs)  by  using  ai  equivalent  value  of  a  which  is  smaller 
than  the  free  stress  value.  This  would  decrease  which  would  increase 
the  dimensionless  drag  paraneter  o^.  As  has  keen  shown,  this  has  the 
effect  of  decreasing  significantly  the  coeputed  velocity  with  which 
the  drops  strike  the  surface.  It  also  appears  to  have  the  effect  of 
increasing  the  surface  density  of  the  drop  phase  but  it  nust  be  renen- 
bered  that 
Pp  =  «1 

where  q  is  the  number  of  drops  per  unit  volune.  If  n  (and  thus  p^)  is 
evaluated  using  the  effective  value  of  a  the  values  of  on  the  wall 
will  decrease  even  though  those  of  increase.  It  should  be  possible 
to  determine  the  effective  values  of  a  appropriate  to  various  flow  con¬ 
ditions  by  matching  numerical  solutions  of  the  type  discussed  in  the 
present  report  with  experimental  results. 

Conclusion 

In  the  present  section  a  numerical  procedure  to  compute  drop-phase 
flow  fields  in  shock  layers  has  been  developed.  A  computer  program  was 
written  to  carry  out  the  calculations.  The  method  was  verified  by  com¬ 
parison  with  an  analytical  solution  for  flow  over  a  thin  wedge.  Numeri¬ 
cal  results  were  computed  and  presented  graphically  for  various  represen¬ 
tative  cases  in  order  to  illustrate  various  parametric  trends.  Further 
results,  corresponding  to  a  particular  experimental  setup  were  found. 
These  will  be  discussed  in  Section  VI. 


(2.46) 


Introduction 


Ibis  part  of  the  research  is  duwctcd  txmrd  the  dnelopnt  of  a 
rational  theory  which  will  adequately  predict  the  relationship  between 
the  inpact  velocity  and  pressure  and  the  crater  size  in  a  given  rain 
awironseat  on  a  given  material-  A  review  of  literature  (for  recent 
comprehensive  surveys,  see  Tlr/nmn  (1967)  and  Cisenberg  (1970))  has 
shorn  that  there  have  been  many  recent  advances  toward  an  understanding 
of  liquid  impact  erosion.  However,  while  an  impressive  amount  of  ex¬ 
perimental  data  has  been  generated,  most  of  the  analytical  models  sug¬ 
gested  have  been  phenomenological  and  not  much  work  has  been  done  toward 
the  development  of  a  physical  theory  for  quantitative  prediction  of 
erosion  damage.  This  is  understandable  in  view  of  the  very  complex 
nature  of  the  damage  mechanism.  The  theoretical  analysis  described 
b&low  is  not  in  itself  expected  to  quantitatively  predict  all  the  rain 
erosion  characteristics  of  materials;  nevertheless,  it  is  hoped  that 
it  will  provide  a  better  understanding  of  the  erosion  mechanism  and 
help  identify  the  parameters  that  are  important . 

The  models  proposed  in  this  study  cue  based  upon  the  assumption 
that  for  high  speeds  of  impact,  the  initial  pressure  developed  is  high 
enough  to  liquefy  the  target  material  in  the  neighborhood  of  the  impact 
area.  Recent  studies  by  Heymann  (1969)  have  shown  that  initial  pres¬ 
sures  of  several  times  the  magnitude  of  "water-hammer  pressure"  are 
possible.  At  sufficiently  high  velocities  of  impact,  therefore,  these 
impact  pressures  are  likely  to  exceed  the  mechanical  strength  of  the 
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target  Material,  making  it  fmfc— i  like  a  liquid  ia  the  i—ffl ate 
neighborhood  of  tke  point  of  impart, 

Sack  a  flea-fliU  iagact  Model  for  ki^ -velocity  impacts  was 
first  coveiicd  by  Cpik  (1936).  Becently  this  Model  was  stwdied  both 
egceiaotally  and  anlytiealty  by  Eagel  (1966,  1967).  Her  analytical 
study,  however,  was  send  empirical,  and  it  involved  an  arbitrary  as¬ 
sumption  with  regard  to  partitiaB  of  tbeesesgy  of  the  inpingjag  drop. 

It  was  ain—ert  arbitrarily  that  half  of  this  energy  is  imparted  to 
the  target  and  the  other  half  is  retained  by  the  drop. 

The  initial  part  of  the  present  research  was  devoted  to  further 
study  of  the  Engel’s  Model  of  rain  erosion.  It  was  hoped  that  it  would 
be  possible  to  use  a  nose  rational  basis  for  the  determination  of 
energy  partition  than  that  used  by  Engel.  It  was  soon  realized,  how¬ 
ever,  that  the  energy  balance  associated  with  liquid-solid  impact  is 
very  complex.  Part  of  the  kinetic  energy  of  the  impinging  drop  goes 
to  the  lateral  outflow,  part  is  dissipated  in  the  shock  waves  passing 
through  the  drop  and  in  the  shearing  associated  with  change  of  liquid 
flow  direction,  and  the  remaining  part  is  dissipated  in  the  target 
material.  Consideration  of  all  these  energy  terns  made  the  problem 
unwieldy,  and  attempts  to  refine  Engel's  Model  were  subsequently 
abandoned . 

Instead  it  was  found  that  incorporation  of  the  concept  of  "appar¬ 
ent  mass,"  used  by  Ludloff  (1967)  in  the  treatment  of  hypervelocity 
impacts  of  solids  on  solids,  into  the  Opik-Engel  model  makes  this  model 
more  amenable  for  application  to  the  rain  erosion  problem.  Furthermore, 
this  concept  allows  the  calculation  of  mass  and  energy  of  the  ejecta 


from  the  crater. 


aritfc  the 


spaa d  jit  a  a  lifd4  wlia  rill  k  eawMnwd.  Ms  ifc!  — i  — - 
taciielj  u.vlyrod  and  In  yielded  good  icndts.  Is  yofetad  oat  tf 
fWrtdr,  Wrtnllf,  ad  Sickel  {1969),  the  cavities  funned  by  the  i»> 
pfnganent  of  high-speed  jets  on  a  liquid  nsriSaoc  bear  a  strikiaf  z m- 
raHaw  to  these  ferae J  cm  solid  aarfaees  by  higfr  wmd  impact. 
dale  high-speed  jets  ace  commmonly  need  to  rcSm  {see, 

for  example,  Hnamitt,  et  al  (1967)),  little  week  has  been  done  to  ex¬ 
ploit  this  analogy. 

Jet  Model 

As  indicated  earlier,  this  aodel  is  based  span  the  assumption 
that  for  ugb  speeds  of  iapaet,  the  initial  pressure  developed  is  high 
enoqgh  to  liquefy  the  target  naterial  allowing  the  iapaet  to  be  cate¬ 
gorized  as  a  fluid- fluid  iapaet.  Based  on  the  arguments  of  Plesset 
and  Chapman  (1971)  and  others,  it  is  assumed,  however,  that  the  high 
contact  pressure  due  to  the  initial  iapaet  decays  very  soon  and  Is  re¬ 
duced  to  the  stagnation  level.  The  calculations  of  Plesset  and  Chapnan 
showed  that  the  duration  of  the  stagnation  pressure  pulse  is  an  order 
of  nagnitude  higher  than  the  period  over  which  the  initial  water  banner 
pressure  acts.  They  further  expressed  the  opinion  that  this  longer 
acting  stagnation  pressure  pulse  nay  be  the  primary  source  of  danage 
caused  by  fluid  impact. 

B  ased  upon  the  above-mentioned  arguments ,  a  simple  idealized  model 
of  steady-state  impingement  of  a  cylindrical  liquid  jet  on  a  liquid 
surface  will  be  used.  The  drop  is  replaced  by  a  cylindrical  jet  in 
order  to  simplify  mathematical  analysis.  This  simplification  is  prompted 
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by  the  fact  tkst  jets  see  oftw  sssd  with  success  to  sisslsts  rail 
otfiiai.  lift  of  its  iopnrfrr  5»  indue tr  i  ■!  pmra»wt  (sack  as 
the  aw  conersioa  proem  is  tk*  steel  iatetxy),  the  jet  inpiug*  •" 

■eat  oa  liquid  surfaces  has  boat  stadied  nteosiwly  both  eqcrieatally 
and  aulytically.  Sane  of  this  work  is  described  below. 

tab  and  Chandrasekhara  (1963)  presented  an  cqeriastal  iawsti- 
gation  of  nonal  penetratioa  of  a  high-«locity  gas  jet  through  a  liquid 
surface,  tab  and  Bawwi  (1965)  extended  this  study  to  a  nee  general 
cue — experimental  study  of  the  iupingenent  of  a  liquid  jet  on  the  surface 
of  a  heavier  liquid.  Sane  300  experinents  were  made  to  measure  the  sizes 
of  cavities  fbrwed  due  to  iapact  of  oil  jets  on  water  and  water  jets  on 
carbon  tetrachloride. 

Frau  the  viewpoint  of  the  present  application,  the  aost  useful  ana¬ 
lytical  study  of  the  problew  is  by  Rossler  and  Stewart  (1968).  They  de¬ 
rived  the  governing  differential  equation  for  the  cavity  profile  using  a 
steady-state  force  analysis.  However,  since  their  eventual  purpose  was 
to  explore  the  instability  of  the  indention  famed  by  the  inpact,  they 
did  not  verify  their  equation  against  experinental  results  such  as  those 
cited  above.  In  what  follows,  a  general  for*  of  the  Rossler- Stewart 
equation,  applicable  to  any  given  pressure  and  shear  distribution  on 
the  liquid  surface,  will  be  derived.  After  testing  the  theory  against 
experimental  results  of  Banks  and  Bhavamai  (1965)  for  the  case  of  iapact 
of  a  liquid  jet  against  a  liquid  surface,  it  will  be  applied  to  determine 
the  crater  sizes  in  the  high-speed  impingements  of  water  and  mercury  drops 
against  metals.  Attempts  will  be  made  to  determine  the  correspondence  of 
surface  tension  with  some  strength  property  of  the  target  material.  Experi¬ 
mental  results  of  Engel  (1959,  1960)  will  be  used  for  this  purpose.  An 


P  =  {WWW  difference  across  tk  interface 
t  =  shear  stress  on  the  interface 
a  =  surface  tension  of  the  target  liquid 
w  =  weight  of  the  inpinging  jet  occupying  the  canity 
B  =  buoyant  force 

For  the  cavity  to  be  in  equilibrium,  the  sue  of  the  vertical  com¬ 
ponents  of  these  forces  oust  vanish.  Hence 
2xrds  p  casfl  ♦  2*Irosin0lp 

=  2xrds  t  sine  i-  2xrdr(pj  -  P2  )gh  +  Jirfrosine]^  =  0 
Here  pj  and  02  are  jet  and  target  densities  respectively,  g  is  the  gravi¬ 
tational  constant,  and  6  is  the  angle  between  ds  and  a  horizontal  line. 

Noting  that  ds  =  dr/cos0  and  taking  the  limit  results  in 

p  -  t  tan0  =  (p2  -  P^Jgh  +  (o/r*Wdr(r  sine)  (3.1) 

Finally,  if  the  trigonometric  relation  tan9  =  -dh/dr  and  a  corresponding 
one  for  sin 8  are  used,  equation  (3.1)  becomes 

h"  +  (hVr)tH-(h')2]  =  (l/o)[(p2-P1)gh-p-Th,][l+T(h’)233/2  (3.2) 

Here  h  and  r  sure  coordinates  describing  the  cavity  profile  (see  Figure  3.1a) 
and  the  primes  denote  differentiation  with  respect  to  r.  Equation  (3.2)  is 
the  governing  equation  for  the  cavity  profile  and  is  to  be  solved  subject  to 
the  following  boundary  conditions: 
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h*  :0itr>0  md  at  p-*  -  bodt  h*  and  h*0 

C^BatuM  (3.2)  is  a  hl^ly  vtfiwl  intar  itcccd  order  diffgtBtial 
eqnadoB.  The  iapvt  infanstioo  needed  for  its  solution  includes  the 
two  densities,  the  surface  tension  of  the  target  fluid,  and  the  pressure 
and  shear  distributions  along  the  cavity  profile.  Unfortunately  no  ana¬ 
lytical  fanes  of  pressure  and  shear  distributions  are  available  for  high¬ 
speed  inpacts.  Approximate  fores  based  on  experinental  result*  will 
therefore  be  needed. 

Before  discussing  the  possible  fanes  of  pressure  and  shear  distribu¬ 
tions  for  the  solution  of  equation  (3.2),  it  will  be  convenient  to  non- 
dinensionalize  the  equation.  The  following  relationships  are  teed  for 
this  purpose: 

h*  =  h/BR,  r*  =  r/BR,  p*/(%>i^),  t*  =  r/dpiV2) 

Here  R  is  the  jet  radius,  V  is  the  impact  velocity,  and  8  is  an  adjust¬ 
able  nondimension al  parameter  chosen  so  that  r*  =  1  will  be  the  location 
of  zero  pressure  point.  The  starred  quantities  are  nondioensional.  The 
dimens ionless  fora  of  equation  (3.2)  can  then  be  written  as  follows : 
h*"  +  (h*/r*)[l  +  (h** )2]  =  [Ah*  -  B(p*  +  i*h*’ )j[l  +  (h*’ )2]3r 
Here  A  =  [(p2-Pi)g82R23/o  and  B  =  (Js p^BRJ/o.  A  and  B  are  both  dimen¬ 
sionless  parameters.  Approximate  forms  of  the  pressure  and  shear  distri¬ 
butions  will  now  be  described  in  terms  of  the  new  dimensionless  variables. 

Pressure  Distribution — In  view  of  the  assumptions  stated  above,  the 
form  of  pressure  distribution  needed  is  that  for  steady-state  impact  of 
a  liquid  jet  on  a  flat  surface.  Banks  and  Chandrasekhara  (1963)  used  a 
normal  distribution 
p*  =  exp(-Kr*2) 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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wbich  served  as  a  fair  approximation  to  the  aeasured  data  of  Gibson  (193%). 
this  distribution,  bowser,  doss  not  take  tbe  possible  negative  pressure 
near  tbe  cavity  edge  into  account.  Rosslen  and  Stewart  (1968)  used  a 
different  distribution  to  fit  Gibson's  data,  a  Modified  font  of  which 
is  given  by 

p*  =  F  sin(l>r*)v/2  (3.7) 

1  ,  0  -  r*  <  1 

where  F  =  <  0.015  (r*-l),  1  -  r*  <  3 

0  ,  3  -  r* 

This  sectionally  continuous  distribution  does  take  the  negative  pressure 
near  the  cavity  edge  into  account.  A  conservation  of  momentum  analysis 
showed  that  the  parameter  6  in  the  definition  of  r*(r*  =  r/0R)  must  be 
equal  to  2.72.  However,  a  value  of  8  =  2.0  provides  a  better  fit  with 
the  experimental  data  of  Gibson.  The  pressure  distributions  given  by 
equation  (3.6)  and  equation  (3.7)  with  8=2.0  and  0  =  2.72  are  compared 
with  Gibson’s  results  in  Figure  3.3a. 

Another  useful  equation  for  pressure  distribution  is  that  by  Leech 
and  Walker  (1966).  They  carried  out  an  extensive  set  of  experiments  and 
fitted  their  results  with  the  following  polynomial: 

p*  =  1  -  3I**2  +  2r*3,  0  =  2.6  (3.8) 

Their  res  Its  however  do  not  differ  significantly  from  those  of  Gibson. 

A  pressure  distribution  which  is  close  to  that  of  Leech  and  Walker  and 
which  takes  negative  pressures  near  the  cavity  edge  into  account  is  given 
by  the  equation 

p*  =  (1  -  r*2)  exp  (-r*2),  0  =  2.6  (3.9) 

The  pressure  distributions  given  by  equations  (3.7),  (3.8),  and  (3.9)  are 
compared  with  each  other  in  Figure  3.3b. 
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Shear  Distribution — Banks  and  Chandrasekhara  (1963)  and  Cheslak, 
et  al  (1969)  concluded  fro*  their  experimental  results  that  the  viscosity 
of  the  inpinging  jet  has  a  negligible  effect  on  stable  cavities  on  liquid 
surfaces.  In  the  impact  of  drops  on  solid  surfaces,  however,  the  shear 
stresses  may  not  be  negligible  because  the  velocity  of  the  lateral  flow 
is  expected  to  be  very  large.  To  take  the  shear  stresses  into  account, 
the  following  approximate  distribution  suggested  by  Rossler  and  Stewart 
(1968)  will  be  used: 

t*  =  (o2 ) ( 0 . 664/Reg ) { 1  -  exp[-(s/R)2]}  (3.10) 

Here  Reg  =  p^aVs/y  is  the  local  Reynolds  number  based  on  the  distance  s 
along  the  cavity  profile,  y  is  the  viscosity  of  the  jet  fluid.  The 
numerical  factor  a  (to  be  called  the  shear  factor)  has  been  introduced 
to  take  care  of  the  possibility  that  the  lateral  flow  velocity  may  be 
larger  than  the  normal  impact  velocity.  In  the  case  of  liquid-solid 
impacts,  several  experimenters  (see,  for  example,  Heymann  (1967))  have 
found  the  lateral  flow  velocity  to  be  several  times  the  velocity  of  impact. 
Numerical  Solution — A  generalized  program  (for  a  listing  of  the  pro- 
j  gram,  see  Hsu  (1972))  usin^  a  fourth-order  Hamming's  method  based  on  an 

|  interval-halving  iterative  technique  was  developed  to  solve  the  govern- 

;  ing  equations  on  a  Xerox  Sigma  6  computer.  This  program  was  used  to 

determine  the  cavity  shapes  and  sizes  in  the  case  of  jet  impingement 
on  liquid  surfaces  as  well  as  liquid  drop  impact  on  solid  surfaces. 

Because  of  the  nonlinear  nature  of  the  governing  equation  and  the 
nature  of  the  boundary  condition  made  convergence  difficult,  several 
subroutines  were  written  to  take  care  of  the  various  difficulties.  The 
final  form  of  the  computer  program  seems  to  provide  fast  and  stable  con¬ 
vergence  in  all  cases  under  consideration. 
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’Uscussion  of  Results — First,  the  effect  of  various  parameters  in 
the  governing  equation  was  investigated.  For  simplicity,  the  center 
line  velocity  of  the  jet  was  assumed  to  remain  unchanged  upon  impact, 
thus  neglecting  the  effect  of  jet  spreading.  The  parameters  used  were 
those  for  impact  of  an  air  jet  on  water.  As  shown  in  Figure  3.4,  the 
cavity  depth  increases  with  increasing  velocity.  When  the  velocity 
exceeds  a  certain  limit,  the  cavity  shape  is  no  longer  shallow,  and  a 
cavity  lip  is  created.  Both  the  cavity  depth  and  the  lip  height  con¬ 
tinue  to  grow  until  a  maximum  cavity  size  is  reached  at  a  critical 
velocity.  Beyond  this  value  the  cavity  depth  starts  to  decrease  while 
the  lip  height  continues  to  increase.  It  might  be  mentioned,  however, 
that  the  pressure  distribution  used  be canes  less  and  less  accurate  as 
the  cavity  gets  deeper. 

In  Figure  3.5,  the  effect  of  increasing  the  shear  factor  a  on  the 
cavity  size  is  presented.  The  shear  factor  seems  to  have  a  more  signif¬ 
icant  influence  on  deep  cavities  than  on  shallow  cavities.  Note  also 
that  the  numerical  solution  is  capable  of  providing  the  cavity  profile 
as  well  as  the  cavity  depth.  Figure  3,5  also  shows  the  presence  of  a 
cavity  lip  in  the  case  of  deep  cavities. 

Next,  the  effect  of  increasing  the  surface  tension  was  investigated. 
As  expected,  the  depth  of  the  cavity  decreases  as  the  surface  tension  is 
increased.  The  effect  of  changing  the  parameter  6  in  the  pressure  dis¬ 
tribution  will  be  discussed  later. 

The  theory  was  then  checked  against  experimental  results  of  Banks 
and  Bhavamai  (1965)  for  impacts  of  oil  jets  against  water  and  water  jets 
against  carbon  tetrachloride.  The  effect  of  turbulent  spreading  of  the 
jet  was  taken  into  account  in  this  case. 
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Figsr«  3.6  presents  the  results  for  the  ease  of  the  Inpact  of  an 
oil  jet  on  water.  The  coordinates  used  are  the  sane  as  those  used  by 
Banks  and  Bhavamai.  The  ordinate  h*^  is  the  cav'ty  depth  ncndimension- 
alized  with  respect  to  the  jet  elevation.  The  abscissa  H*  represents 
nondinensional  jet  nonentus.  The  experimental  points  of  Banks  and 
Bhavanai  for  various  nozzle  elevations  are  plotted  on  this  graph.  Also 
plotted  are  the  calculated  results  of  the  present  theory  using  pressure 
distribution  equation  (3.7)  with  B  =  2.72  and  B  =  2.0.  In  vies  of  the 
considerable  scatter  in  the  experimental  data,  the  agreement  seems  to 
be  satisfactory. 

An  empirical  equation  obtained  by  Banks  and  Bhavamai  which  fits 
their  experimental  data  is  compared  with  the  present  theory  in  Figures 
3.7  and  3.8  for  oil-water  and  water-carbon  tetrachloride  impacts  respec¬ 
tively.  The  ordinate  used  here  represents  a  normalized  cavity  depth 
defined  in  Banks  and  Bhavamai  (1965)  or  Hsu  (1972).  Specific  compari¬ 
sons  with  experimental  results  showed  that  the  higher  value  of  3  is 
better  for  cases  of  deep  penetration  while  lower  values  of  3  are  more 
suitable  for  shallow  cavities. 

This  model  was  also  applied  to  the  case  of  the  impact  of  a  water 
drop  on  water  surface.  This  case  was  considered  by  Engel  (1966)  as  a 
model  for  high-speed  liquid-solid  impingements.  In  Figure  3.9,  her 
experimental  ■results  are  compared  with  those  computed  from  the  present 
theory.  The  agreement  is  satisfactory. 

The  fluid-fluid  impact  model  was  then  applied  to  the  impact  of 
liquid  drops  on  solid  surfaces.  The  premise  on  which  this  model  is 
based  has  already  been  discussed  earlier.  In  this  analysis,  the  pres¬ 
sure  distributions  given  by  both  equations  (3.7)  and  (3.9)  were  used, 
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but  no  significant  difference  was  noted.  As  far  as  the  shear  distribution 
given  by  equation  (3.10)  is  concerned,  it  was  found  that  increasing  the 
shear  factor  a  to  a  value  up  to  4  does  not  significantly  change  the  re¬ 
sults;  therefore,  a  was  fixed  at  a  value  of  1.  In  fact,  it  was  found  that 
the  shear  term  in  equation  (3.5)  is  small  compared  with  other  terms  in 
the  equation.  Several  workers  (see,  for  example.  Smith  and  Fyall  (1969)) 
have  found  from  their  experiments  that  the  shearing  action  resulting  from 
the  high-speed  radial  flow  is  not  very  significant. 

One  of  the  difficulties  in  applying  the  proposed  model  to  liquid- 
solid  impacts  is  that  of  relating  the  surface  tension  parameter  a  to  some 
strength  property  of  the  solid.  It  was  fortunately  found  that  the  use  of 
a  simple  relationship 

S  =  o/R,  (3.11) 

where  S  is  the  dynamic  yield  strength  of  the  target  material  and  R  is  the 
drop  radius,  gives  good  correspondence  with  the  available  experimental 
results . 

Figures  3.10,  3.11,  3.12,  and  3.13  present  the  calculated  results 
for  the  impact  of  2-mm  water  and  mercury  drops  on  copper  and  aluminum 
targets  at  high  velocities.  Engel  (1960)  carried  out  experiments  to 
measure  pit  depths  in  such  impacts.  Mercury  drops  were  used  because 
they  are  capable  of  imparting  large  momenta  due  to  their  high  densities. 
Experimental  results  of  Engel  are  also  plotted  on  the  figures  for  com¬ 
parison.  Calculated  results  seem  to  be  surprisingly  good  in  view  of 
the  crude  nature  of  the  model  used.  Because  of  the  assumption  that  the 
target  material  behaves  as  a  fluid,  the  model  should  be  more  applicable 
at  higher  values  of  the  momentum.  This  is  apparent  from  the  figures; 
the  correspondence  between  calculated  and  measured  pit  depths  is  better 
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for  the  —retry  impact  than  fop  the  water  impact.  Also,  la  the  case  of s 

:  *  f 

water  Iapact,  the  calculated  results  get  closer  to  the  seasoned  data  as 

S  i 

»  s 

tiie  iaplng— eat  velocity  increases.  _ 

Approximate  Solutiao--Examination  of  the  relative  importance  of  the 
*  1 
various  ‘pan— stars  in  the  nondinensicnal  fore  of  the  governing  equation 

l 

(equation  (3.5))  revealed  the  .possibility  of  a  siaple  approximate  solu¬ 
tion.  In  the  case  of  liquid-solid  inpacts  a  is  very  large  and  conse- 

.  I  j 

quently  the  f  nondinensicnal  par— etir  A  is  very  snail.  The  first  tern 

1  i 

on  the  right-hand  side  is  therefore  negligible  in  comparison  with  otter 

terns.  As  suggested  earlier,  the  shear  tern  makes  an  insignificant  con- 

«  » 

tribution.  In  addition,  the  cavities  are  expected  to  be  shallow,  there*' 
fore  (h*')2  '«  1.  With  these  approximations ,  equation  (3.5)  simplifies 
to 

■  !  I 

h*"  +  h*'/r*  =  -Bp* 

An  exact  solution  to  equation  (3.12)  can  be  obtained  if  the  approxi- 

I 

mate  form  of  the  pressure  distribution  given  by  equation  (3.:9)  is  used. 

i 

The  solution  of  equation  (3.12)  subject  to  the  proper  boundary  conditions 
(equation  (3.3))  is  then,  simply,  . 

. #/  l 

h*  =  (B/4)  exp  (-r*2)  ‘  1 

*  l 

The  maximum  cavity  depth  which  occurs  at  r*  =  0  is  given,  in  dimensional 

’  .  '  •  ,  i 

i  ’ 

form,  by  the  following  expression: 

W  *  Pl*W/8o 

The  approximate  solution  given  by  equation  (3.i4)  is  compared  with 
the  numerical  solution!  of  the  entire  equation  and  with  the  experimental 

results  of  Engel  in  Figures  3,14  and  3.15  for  impacts  of  mercury  drops 

‘  \ 

against  2024-0  aluminum  and  steel  respectively.  It  can  be  concluded  that 
the  approximate  solution  is  quite  adequate.  Similar  correspondence  was 
found  in  other  cases  not  presented  here. 


(3.12) 


(3.13) 


i 

(3.14) 


I 


-36- 


It  nay  be  coo  eluded  fra  the  above  results  that  the  jet  Model 
described  above  is  capable  of  providing  good  results,  at  least  in  the 
case  of  liquid  inpact  on  Metals.  Unfortunately ,  experimental  results 
in  a  fom  suitable  for  checking  the  validity  of  the  Model  were  not 
available  for  other  Materials,  such  as  cerand.cs.  However,  the  present 
results  are  significant  in  two  respects.  First,  the  results  may  be 
directly  applicable  because  sane  designs  of  radones  suggest  that  Metals 
nay  be  incorporated  iu  the  attachment  to  the  afterbody  of  the  vehicle 
or  in  the  nose  as  a  protection  against  rain  erosion.  Microwave  engi¬ 
neers  have  shown  that  such  a  use  of  Metals  is  possible  without  serious 
detrinent  to  radar  performance.  Second,  the  development  of  a  reliable 
theory  mist  be  based  upon  materials  on  which  the  test  data  is  repeatedly 
reproducible.  Metals  are  such  Materials. 

Apparent  Mass  Model 

The  apparent  mass  model,  just  like  the  jet  model,  is  a  rough, 
approximate  hydrodynamic  model.  This  model,  based  on  a  concept  first 
used  by  Ludloff  (1967)  in  connection  with  his  study  of  hypervelocity 
impact,  is  characterized  by  the  following  assumptions: 

1.  Upon  impact  a  fraction  of  the  initial  kinetic  energy  is 

expended  in  liquefying  the  surrounding  target  material.  The 
remainder  of  the  initial  kinetic  energy  causes  impulsive 
motion  in  the  liquefied  mass  such  as  that  which  occurs  when 
a  body  impacts  on  the  free  surface  of  a  liquid.  Such  a 
motion  can  truly  be  characterized  as  irrotational  (see 
Batchelor  (1967)).  The  theory  of  flows  with  free  surfaces 
suggests  that  the  concept  of  "apparent"  mass  may  be  useful  in 
describing  the  flow  of  the  target  material. 
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2.  Since  tbe  drop  is  expected  to  be  flattened  by  iapact,  tbe 
velocity  field  of  tbe  target  material  is  similar  to  that 
around  a  disc  translating  tiragb  a  fluid. 

3.  Tbe  arterial  resists  the  fowaation  of  the  crater  by  a  force 
equal  to  tbe  prodt  ,t  of  sone  strength  property  S  of  the  eater- 
ial  and  the  surface  area  of  the  crater. 

Wien  a  body  of  nass  N  enters  a  liquid,  there  is  a  sudden  reduction 
of  its  velocity.  This  is  ascribed  to  an  apparent  addition  to  its  nass 
by  tbe  nass  of  the  liquid  set  instantaneously  in  notion.  This  additional 
nass  is  called  the  apparent  nass. 

In  the  present  case  it  will  be  ass  used,  for  simplicity,  that  the 
crater  is  henispherical  and  renains  in  that  shape  as  it  grows  with 
increasing  radius  r.  The  equation  of  notion  for  the  target  fluid  can 
then  be  written  as 
d/dt[(M  +  M* )r]  +  2*r*S  =  0 

Here  H  is  the  mass  of  the  drop  and  M'  is  tbe  apparent  mass.  For  the 
motion  of  a  disc,  the  apparent  nass  is  given  by  (see  Batchelor  (1967)) 

M‘  =  (8/3)p2c3 

where  p2  is  the  density  of  the  target  material  and  c  is  the  radius  of  the 
disc. 

If  M  and  H*  are  constant,  equation  (3.15)  can  be  integrated  over  r 
to  obtain  the  energy  equation 
(l/2)(M  +  M»  )r*  +  (2/3)irSr3  =  (1/2)MV2 

Here  the  condition  r  =  V  when  r  =  0  has  been  applied.  V  is  the  impact 
velocity.  The  maximum  cavity  depth  is  the  value  of  r  when  r  =  0  and  is 
given  by 

fynax  =  (3MV2/4rS)1/3 


(3.15> 


(3.16) 


(3.17) 


(3.18) 


follow: 


E  =  (3.19) 

Here  t _ is  the  time  required  to  achieve  a  cavity  depth  of  aapitode 

h_1 _ The  time  can  easily  be  obtained  by  integrating  equation 

(3.17). 

The  solution  outlined  above  was  tested  against  known  experimental 
data.  Although  the  moaex  seew  to  be  qualitatively  correct,  quantitative 
correspondence  was  lacking.  This  is  hardly  surprising  in  view  of  the 
very  rough  nature  of  the  aodel. 

It  seems  that  the  model  may  be  improved  by  two  additional  considera¬ 
tions.  First,  no  provision  has  been  made  in  the  model  for  ejected  material 
during  crater  formation .  This  can  be  done,  as  suggested  by  Ludloff  (1967), 
by  specifying  an  additional  integral  condition  based  upon  experimental 
results  about  ejecta.  Second,  the  expression  far  the  apparent  mass  (equa¬ 
tion  (3.16))  should  be  modified  to  take  into  account  the  growth  of  the 
crater.  The  following  expression  is  suggested: 

M'  =  ep2r3  (3.20) 

where  e  is  an  adjustable  parameter. 

Equation  (3.15)  can  then  be  written  as 

(d/dt)[ (H  +  ep2r3)r]  +  2*r2S  =  0  (3.21) 

This  equation  has  to  be  solved  numerically.  Efforts  are  currently  under¬ 
way  to  determine  the  solution  and  to  test  it  against  experimental  results. 

The  apparent  mass  model  deserves  further  consideration,  because,  as 
opposed  to  the  jet  model,  it  considers  the  process  to  be  time-dependent. 


High  up— d  calliriwi  ritk  objmcts 


defans  Is  space  or  dettf  K«i ,  asd  nladfops  Is  tfa  4lu<{)hut  pose 
potential  haaurds  to  spacecraft,  aianiles,  or  airplanes. 

A  hi^Helodty  impact  gene  rally  creates  a  enter  is  the  structure, 
driving  a  strong  shock  wave  into  it.  If  the  impacted  structure,  or 
"target,"  is  sufficiently  this,  a  picture  will  result.  If  tie  target 
is  relatively  thick,  the  shock  will  rapidly  decay  into  an  elastic  stress 
wave.  Hben  such  a  wave  encounters  a  free  surface  it  is  reflected,  gen¬ 
erally  as  a  tensile  wave,  and  its  anplitude  way  be  of  sufficient  Magni¬ 
tude  to  produce  fractures  near  the  rear  surface.  Such  fractures  may 
appear  as  granular  cracks  near  the  surface,  as  rear  surface  bulges,  or 
as  a  ooMplete  detachment  of  target  Material,  creating  a  shrapnel  effect. 

An  example  of  each  is  shown  in  Figure  4.1  far  copper,  altwunun,  and 
steel  targets. 

The  fractures  in  a  transparent  acrylic  (methyl  methacrylate)  resin 
are  also  shown  in  Figure  4.1.  Although  this  material  (trade  names  of 
Lucite  and  Plexiglas)  would  probably  never  be  used  for  the  actual  struc¬ 
ture,  it  is  excellent  for  the  study  of  stress  waves  and  fracture  as  both  the 
waves  and  the  fracturing  process  can  be  photographed.  These  can  then  be 
compared  with  analytical  conclusions  (Kins low  (1964)). 

Stress  Relations  for  Spherical  Waves 

Spherical  dilataticnal  wave  propagation  in  a  homogeneous ,  iso.Topic 
material  can  be  specified  _>y  the  equation 
32$/3t2  =  (^[a^/Jr2  +  (2/r>3#/3*9 


(4.1) 
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Mhoc  •  is  a  scalar  diiplac— art  potutlal,  is  the  save  atlodtr,  aai 
t  is  tiae.  Particle  ditplacr— nt  (a)  mi  velocity  (»)  as*  specified  by 
the  relatioas 

a  =  as/ar  and  f  *  »/X  (4.2) 

dare  r  denotes  tbe  radios  vector  brae  the  point  of  projectile  impact. 

The  radial  and  tafatial  stresses  are  given  by  the  relations 

oT  *  (1  ♦  2e)(3B/ar)  ♦  2>(u/r)  (4.3) 


ot  s  ot  =  l(3d/jr)  t  2C1  ♦  y)(u/r)  (4.4) 

1  2 

where  1  and  y  are  the  lad  constants  and  are  related  to  Young's  Modulus 
(E)  and  Poisson's  ratio  (v)  as  follows: 

X  =  vE/Cl  *  v)(l  -  2v)  (4.5) 

V  =  E/2(l  ♦  v)  (4.6) 


Model  for  Generating  Stress  Waves 

The  nathenatical  nodel  for  generating  spherical  elastic  waves  is 
that  described  by  tbe  author  in  the  reference.  Kins  low  (1963).  To  suw- 
aarize  briefly,  it  is  assuned  that  there  is  a  hollow,  hemispherical 
cavity  in  the  target  with  its  center  at  the  point  of  impact  and  that  a 
tine-varying  pressure  or  forcing  function  is  applied  to  this  imaginary 
cavity  surface,  generating  stress  waves  in  the  material.  The  pressure 
applied  to  the  surface  of  ’he  cavicy,  having  a  radius  ro,  is  an  impulse 
described  by  the  relation 
p/p0  =  K(e  °lt:  -  e  °2t) 

where  pQ  is  the  maximum  pressure,  and  Oj  are  decay  constants,  t  is 
elapsed  time,  and  K  is  a  constant.  By  the  proper  choice  of  values  of 
and  o2  various  wave  forms  can  be  generated.  The  values  for  describ¬ 
ing  the  forcing  function  may  be  selected  to  approximate  liquid  impact. 


(4.7) 


ical  cavity,  a  solution  of  tin 


*  s  i  iKiVoT0}/{pT<ti  *  iao~  ia)2J)(  ~*iace 

i=l 

♦  e~*°T  (  1  *  ((oo-im)/^))2)1*  cos(wor  -  tan_1(«0-i»)/(w0)))) 

(*.0) 

where 

r  *  t  -  ((r-r0)/(Cj)) 

(*.9) 

*o  =  <Cl/ro)((1"2w)/(1"v)) 

{*.10) 

wG  =  (C1/r0H((l-2v)^)/(l-v)) 

(*».n) 

Reflected  Stress  Haves 

The  simplest  example  of  the  reflection  of  an  elastic  dilatation 
wave  fron  a  free  surface  occurs  when  the  wave  strikes  norwl  to  the 
surface.  Since  the  resulting  surface  stress  oust  be  zero,  a  compression 
wave  oust,  therefore,  always  be  reflected  as  a  tensile  wave,  and  a  tensile 
wave  aust  be  reflected  as  a  conpression  wave.  When  such  a  wave  strikes  a 
free  surface  at  an  oblique  angle,  the  situation  is  nuch  no  re  conplex.  Mot 
only  will  waves  of  dilatation  be  reflected  but  there  will  also  be  generated 
distortional  or  shear  waves.  From  an  ai.ilysis  of  these  wavas  reflected 
from  a  plane  surface  (Kin slow  (1969)), it  seems  that  only  the  incident  and 
reflected  dilatation  waves  need  to  be  considered  in  determining  the  maximum 
tensile  stress  except  at  points  some  distance  from  the  normal  axis  where 
fractures  are  not  likely  to  occur.  For  twis  reason,  the  amplitude  of  the 
reflected  shear  waves  are  not  computed  in  this  preliminary  analysis ,  al¬ 
though  their  existence  is  recognized  and  taken  into  account  in  the  compu¬ 
tation  of  the  magnitude  of  the  reflected  dilatation  waves. 


Relative  values  of  mimtuT  stmui  soon  rftw  tte  wave 

reflection  fna  a  plate  surface  are  shown  in  Figure  4.2.  Coapwtiie 
stresses  are  considered  as  positive  and  tensile  stresses  as  negative. 

At  the  tine  shorn,  the  nasi  nun  tensile  stress  has  a  value  of  -61  units. 
It  can  be  seen,  however,  that  this  is  not  the  sasiana  tension  that  will 
be  created  as  at  a  slightly  later  tine  this  reflected  tensile  wave  will 
cnabine  with  the  tensile  "tail"  of  the  incident  wave  to  generate  a  each 
higher  tensile  stress  at  a  somewhat  greater  distance  from  the  rear  sur¬ 
face. 


Stresses  Resulting  iron  Multiple  Impacts 

The  coordinate  systen  is  given  in  Figure  4.3.  The  locations  of 
inpacts  on  a  plane  surface  (the  x-y  plane)  are  designated  as  — 

and  are  given  by  the  coordinates  (*on»y0n).  ^  location  of  any  point  P  at 

which  the  principal  stresses  are  desired  is  specified  by  the  coordinates 
(x,y,z).  For  the  computation  of  stresses  resulting  fron  Iwpact  at  Cn,  the 
coordinates  are  enployed,  where 

*n  =  x  -  xon 

yn  =  y~  yon 
^  = z 

The  location  of  P  may  also  be  specified  by  the  coordinates  (rn,0n,#n). 

The  following  derivation  of  the  components  of  stress  at  point  P 
are  relative  to  the  0n  coordinate  system.  For  the  sake  of  simplicity 
the  subscript  n  is  omitted  from  the  derivations  and  in  the  six  cartesian 
components  of  stress,  the  symbol [  3n  indicates  that  all  values  enclosed 
by  the  bracket  are  relative  to  that  system. 

The  relations  between  the  (x,y,z)  and  the  (r,6*$)  coordinates  are 


-) 


(4.12) 

(4.13) 

(4.14) 


an. 


x  -  r  slot  cos# 
y  -  r  sin*  sine 
z  =  r  cos# 

The  direction  cosines  of  the  positive  r-axis  bate  components  in  the 
(x,y,z)  systee  as  follows 

cos(x,r)  =  3x/3r  =  sin#  cosd 
cos(y,r)  =  Jy/Jr  =  sin#  sinfi 
cos(z,r)  s  3z/ar  =  cos# 

Sinilarly,  the  direction  cosines  of  -die  t0  -  gt^-axis  are 
respectively 

cos(x,t0)  *  (1/r  sin8)(3x/36)  =  -sine 
cos(y,t0)  =  (1/r  sin0)(3y/38)  =  cos0 
cos(z,t0)  =  0 

and 

cos(x,t0)  s  (l/r)(3x/3#)  =  cos#  cos8 
cos(y,t0)  =  (l/r)(3y/3#)  =  cos#  sine 
cos(z,t0)  =  (l/r)(3z/3#)  =  -sin# 

where 

r  =  /x*  +  y^  +  zs 
sin#  =  /x2  +  y2/r 
cos*  =  z/r 
sinfl  =  y/Zx2  +  y2 
cose  =  x/Zx2  +  y2 

As  the  spherical  wave  propagates  through  the  material,  the  radial 
and  tangential  stresses  having  values  given  by  equations  (4.3)  and  (4.4), 
respectively,  are  shown  acting  on  an  elemental  volume  in  Figure  4.4.  At 


(4.16)  !,1’ 
/ r  * 

(4.17)  j 


(4.18) 

(4.19) 

(4.20) 

(4.21)  | 

i 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

.(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 
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any  point,  P,  and  for  any  tine.  The  six  cartesian  components  of  stress 
relative  to  the  [o,x,y,z]n  coordinate  system  may  be  computed  from  the 
following  relations. 

s  [ot0cos2(x,tQ)  +  ati cos2(x,t^)  +  orcos2(x,r)3n 
Co^  *  Cot4>cos2(y,t^)  +  ar cos2(y,r)  +  ot0cos2(y,t0)3n 
Coz]n  *  Corcos2(z,r)  +  at0cos2(z,t0)  +  ot$cos2(z,t^)3n 
tTxy3n  =  C°t0cos(x,t0>cos(y,te)  +  ot$cos(x,t$)cos(y,t$) 

+  arcos(x,r)cos(y,r)]n 

Ct  L  =  [o  .cos(y,t. )cos(z,t.)  +  o  cos(y,r)cos(z,r) 
yg  n  tv  ♦  #  r 

+  ot0cos(y,t0)cos(z,ta)]n 

[Tzx3n  =  torcos(z,r)cos(x,r)  +  at(J(cos(z,t4))cos(x,t(J)) 

+  ot0cos(z,te)cos(x,t0)]n 

These  components  of  stress  resulting  from  each  of  the  several 


impacts  can  now  be  added  as  follows: 


°x  =  Cox]l  +  C°x]2  +  +  *  •  ' 

°y  =  Layh  +  Cay]2  +  C°y]3  *  ‘  *  * 

°z  =  t°z]l  +  [oz32  +  ^°z^3  +  *  *  * 

Txy  "  ^Txy^l  +  ^Txy^2  +  ^Txy^3  +  •  •  • 

x  =  [t  ]  +  [t  3_  +  [t  3  +  .  .  . 
xz  xz  1  xz  2  xz  3 

T„_  =  [r  3  t  [r  ],  +  [x  3.  +  .  .  . 

yz  yz  1  yz  2  yz  3 

These  six  components  of  stress  are  now  substituted  in  the  following 


(4.32) 

(4.33) 

(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 

(4.40) 

(4.41) 

(4.42) 

(4.43) 


(4.44) 
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Stresses  Resulting  from  Two  Liquid  Impacts 

toe  of  the  problems  that  must  be  dealt  with  in  the  analysis  of 
radorae  damage  resulting  from  rain  impact  is  the  effects  of  multiple 
drops  and  the  interaction  of  the  stresses  developed  by  the  individual 
drops. 

When  a  rain  field  is  described,  the  probable  spatial  and  time 
distribution  of  drops  interacting  with  a  radome  surface  moving  at  a 
specified  velocity  can  be  determined.  The  effects  of  the  distances 
between  two  drops  and  the  times  between  impacts  are  analyzed. 

The  stresses  computed  are  relative  values  only  and  dimensionless 
units  of  time  and  dimensions  are  employed.  The  radome  thickness  is 
denoted  by  T,  the  stress  by  <T,  and  the  time  by  T.  A  pressure  wave 
typical  of  high-velocity  impact  is  employed. 

First,  two  impacts  at  the  same  time  are  considered.  The  maximum 
tensile  stress  resulting  from  two  impacts  will  be  generated  either 
directly  beneath  the  points  of  impact  or  on  an  axis  midway  between 
these  points.  Figure  4.5  gives  the  values  of  the  tensile  stress  for 
different  distances  between  the  impact  points.  The  maximum  tensile 
stress  for  a  single  impact  has  a  relative  value  of  cT  =  115.  If  two 
drops  are  very  close  together  at  the  same  time,  the  resulting  stress 
will  obviously  be  twice  this  value,  or  230.  As  the  distance  between 
drops  is  increased,  the  developed  stress  will  decrease.  In  this  case 
of  simultaneous  impacts,  it  can  be  seen  that  the  maximum  stress  occurs 
along  the  axis  between  the  two  points.  At  a  distance  of  about  0.8  T, 
the  stresses  directly  beneath  the  impact  points  have  decreased  to  about 
that  resulting  from  a  single  drop,  but  the  stress  on  the  z-axis  has  a 
value  of  almost  200.  The  stress  on  this  axis  continues  to  decrease  as 
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the  distance  between  drop6  increases.  At  a  distance  of  approximately 
1.8  T,  the  stress  on  this  axis  has  decreased  to  the  value  created  by  a 
single  inpact.  This  neans  that  if  the  distance  between  two  drops  is 
greater  than  twice  the  thickness  of  the  target,  the  maximum  stress  can 
probably  be  determined  by  computing  the  stress  developed  by  a  single 

drop  and  that  the  interaction  of  the  stress  waves  can  be  neglected. 

. 

The  effect  of  tine  (At)  between  drops  will  new  be  considered.  The 
stresses  developed  by  two  drops  impacting  the  same  point  are  shown  in 
figure  4.6(a).  As  pointed  out,  the  tensile  stress  will  have  a  value 
of  '30  for  zero  tine  between  impacts.  As  tve  time  increases,  the  stress 
decreases  up  to  a  relative  time  of  0.50,  beyond  which  the  stress  is  the 
sane  as  for  the  impact  of  a  single  drop.  If  the  distance  between  drops 

I 

is  0.4  T,  the  tensile  stress  for  various  times  (At)  is  given  in  Figure  1 

4.6(b).  The  first  impact  is  and  the  later  impact  is  P^.  The  other 
curves  give  this  information  for  various  distances  between  impacts. 

It  should  be  noted  that  these  results  are  based  on  an  assumed 
material  and  for  a  typical  stress  wave.  The  results  would  not  be  the 
same  for  different  materials  and  impact  conditions. 

Stresses  Resulting  from  Single- Drop  Impact 

When  a  collision  occurs  between  a  fast-moving  surface  and  a  water- 

I 

drop,  a  high  pressure  is  first  generated  at  the  point  of  contact.  This  j 

is  followed  by  a  lateral  flow  of  water  at  very  high  velocity.  The  pres-  ! 

sure  distribution  at  the  liquid-solid  interface  and  the  times  and  dura¬ 
tions  of  the  pressures  at  various  points  on  the  surface  have  not  yet 
been  determined  except  for  rather  low  impact  velocities  such  as  that 

being  done  by  Huang  (1971-72).  Similar  data  for  high  Mach  numbers  is  j 

\ 
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perbaps  the  greatest  need  in  the  analysis  of  this  problem  at  the  present 
time.  Both  analytical  and  experimented  methods  are  being  employed  in 
attempts  to  obtain  this  information.  For  the  relatively  low- velocity 
impacts,  Huang's  computations  indicate  that  the  pressure  is  a  maximum 
at  the  center  of  the  impact  area,  decreasing  as  the  distance  from  the 
center  increases.  Some,  including  Herbert  (1970),  believe  that  the 
pressure  is  not  a  maximum  at  the  center,  but  that  the  maximum  pressure 
is  near  the  periphery  of  the  drop.  Others  believe  that  cavitation  occurs 
because  of  the  high  jet  velocity  creating  tensile  stresses  at  the  surface. 

Due  to  the  non-uniform  time-varying  pressure  distribution  generated 
by  liquid-drop  impact,  it  is  not  possible  to  accurately  determine  the 
stress  wave  characteristics  by  a  single  stress  wave  source.  The  method 
of  multiple  impacts  derived  in  a  previous  paragraph  will  be  used  by  con¬ 
sidering  a  cluster  of  wave  sources  as  shown  in  Figure  4.7  for  each  drop 
impact.  The  forcing  function  to  be  applied  to  each  imaginary  cavity  will 
be  determined  by  the  pressure- time  history  over  the  target  surface.  This 
will  result  in  a  large  number  of  wave  fronts  which,  according  to  huygens* 
principle,  will  combine  to  form  a  single  wave.  The  stress  distribution 
created  by  this  wave  can  be  computed  for  various  times.  The  maximum 
pressure,  tension,  and  shear  stresses,  as  well  as  the  displacements, 
strains,  and  strain  rates  can  be  determined. 

The  locations  of  fractures  in  a  Perspex  plate  resulting  from  the 
impact  of  a  water  jet  is  shown  in  Figure  4.8(a).  This  was  taken  from  one 
of  the  many  excellent  photographs  by  Brunton  (1967)  of  Cambridge  Univer¬ 
sity.  It  shows  a  ring  crater  on  the  surface,  a  fracture  directly  under 
the  center  of  the  impact  area,  and  a  large  fracture  near  the  rear  surface. 
Not  taowing  the  actual  pressure  distribution,  various  values  were  assumed 


43* 
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<■  =  ,  .  ,  « 

and  the  maximum  tensile  stresses  were  computed  throughout  .the  target. 

The  areas  where  these  stresses  exceeded  a  specified  value  are  shown  in 
Figures  4.8(b)  through  (g).  '  The  similarity  to  the  observed  fractures  , 

•’  1  ,  *  I 

may  be  seen.  The  stresses  not  only  depend  upon  the  maximum  pressure 
distribution  but  the  time  rate  of  the  spread  of  these  pressures  from 
the  point  of  first  contact.  This  is  illustrated  by  Figures  4.8(e)  and 

I  *  *  *  * 

(f).  The  same  maximum  pressure  distribution  was  assumed  in  both  in-  • 

*  t  i 

stances  and  only  the  times  were  different.  These  shaded  areas  would 

l  .  ! 

represent  fractures  only  if  a  critical  stress  criterion  of  fracture 

i 

applies . 

%  assuming  various  combinations  of  pressures  and  times* it  may 
be  possible  to  obtain  a  numerical  solution  thai  will  match  the  experi¬ 
mental  results.  If  this  is  accomplished*  the  pressure-time  values 

'  » 

used  in  the  computation  should  represent  the  true  values  for  that 

, 

particular  impact.  , 

This  method  of  multiple  stress  wave  sources  can  also  be  applied 

’  i 

for  the  case  of  oblique  impact  if  the  pressure-time  variations  are 

.1  ,1 


known. 


! 
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SECTION  V. 

PRELIMINARY  EXPERIMENTS 

Liquid  Jet  Accelerator 

A  device  similar  to  the  one  constructed  at  the  Cavendish  Laboratory, 

University  of  Cambridge,  and  described  by  Fyall  in  Chapter  8  of  the 
Radome  Engineering  Handbook,  edited  by  Walton  (1970),  has  been  constructed 
and  used  to  accelerate  small  amounts  of  water  to  a  relatively  high  veloc- 
ity.  The  accelerator  is  shown  in  Figure  5.1.  The  nylon  piston  is  in¬ 
serted  in  the  chamber  to  almost  the  aid;  a  few  drops  of  water  are  inserted 
in  the  orifice  with  a  hypodermic  needle;  and  the  piston  is  accelerated  by 
means  of  a  0.22-caliber  rifle  bullet  striking  its  end.  For  the  dimensions 
shown,  a  bullet  velocity  of  750  ft/sec  accelerates  the  water  to  a  velocity 
of  2200  ft/sec.  Photographs  of  the  impact  range  (constructed  in  a  surplus 
iron  lung)  are  shown  in  Figure  5.2.  The  jet  coining  from  the  nozzle  is 
shown  in  Figure  5.3.  Variations  in  the  piston  and  nozzle  design,  the 
amount  of  water,  and  the  bullet  velocities  will  be  made  in  attempts  to 
reach  higher  velocities  and  to  more  accurately  simulate  raindrops. 

i 

Impact  Dangge  to  Lucite  Targets 

I 

Figure  5.4  shows  the  ring  fractures  on  the  surface  of  a  0.25-inch  1 

Lucite  target.  The  center  portion  of  the  surface  was  not  affected.  The 
lower  photograph  on  the  first  page  is  of  a  sectioned  target  so  that  the 
depth  of  fracture  could  be  seen.  Ring  craters  resulting  from  lower 
velocities  are  shown  on  the  second  page  of  the  figure.  These  ring 
craters  are  similar  to  the  ones  photographed  by  Brunton  (1967)  and  by 
Fyall  (1969).  Side  views  of  targets  are  shown  in  Figure  5.5  in  which 
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the  near  surface  fractures  can  be  seen,  there  were  also  fractures  below 
the  center  of  the  point  •'f  inpact  but  are  barely  distinguishable  in 
these  photographs. 

A  vertical  silling  machine  and  a  dial  gage  were  used  as  a  profilaueter 
for  seasoning  crater  depths.  Figure  5.6  shows  the  craters  on  a  noseoone 
being  Measured,  the  photograph,  cross-section,  and  depth  contours  of  a 
crater  in  bucite  are  shown  on  Figure  5.7. 

Pasage  to  Slip-Cast  Fused  Silica  Targets 

The  front  and  rear  surfaces  of  the  SCFS  targets  were  painted  with 
blade  ink  before  inpact  so  that  the  extent  of  any  fractures  could  be 
easily  seen.  The  photographs  at  the  top  of  Figure  5.6  show  the  damage 
to  a  0.25-inch  target.  It  was  completely  penetrated.  The  fracture  as 
viewed  from  the  back  is  interesting.  The  lower  photographs  of  this 
figure  show  the  damage  to  a  slightly  thicker  target.  The  familiar  ring 
crater  is  seen  on  the  front  and  a  crack  is  detected  cn  the  real'.  Cross- 
sections  of  this  crater  are  shown  in  Figure  5.9.  The  maximum  depth  is 
about  0.0085  inch.  Unlike  the  ring  crater  in  Lucite,  the  center,  appar¬ 
ently  undamaged  portion,  shows  a  permanent  deformation.  It  may  be,  how¬ 
ever,  that  the  material  is  not  compressed  but  that  the  crack  whose 
boundaries  are  seen  on  the  rear  surface  extends  internally  to  the  front 
crater  and  that  this  entire  cone  of  material  is  pushed  backward.  Crater 
contours  are  shown  in  Figure  5.10. 

It  appeared  that  the  fracture  of  the  0.40-inch  target  shown  in 
Figure  5.11  was  very  similar  to  the  last  one  discussed  except  for  the 
rear  surface  spall  directly  behind  the  inpact  point  and  that  the  frac¬ 
ture  approximately  1  inch  from  the  center  was  somewhat  more  pronounced. 
Shortly  after  the  first  photographs  were  made,  the  entire  rear  portion 
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shwD  is  tiie  lest  two  photographs  of  this  fifSN  became  detached.  Qpoo 
an  nr—inatioo  of  these  fractures,  Figure  5.12  was  drown.  The  photo¬ 
graphs  of  Figure  5.13  show  an  edge  view  of  the  target  just  before  and 
following  impact.  The  spall  aaterial  can  be  seen.  Surface  rtmiigr 
resulting  from  oblique  inpacts  is  shewn  in  Figure  5.1*. 

A  Possible  Explanation  of  the  Target  Fracture 

The  fractures  indicated  by  A,  B,  aid  C  of  Figure  5.12  are  of 
interest  as  either  of  then  would  probably  cause  greater  daaage  than 
does  the  ring  crater.  The  crack  C  is  of  special  importance  for  the 
reason  that  it  creates  the  largest  spall,  and  hence  wore  structural 
weakening,  than  either  of  the  others  would  have  produced. 

Initially  the  path  of  crack  A  is  determined  by  the  tensile  stresses 
around  the  i Mediate  impact  area.  As  suggested  by  Yoffee  (1951)  the 
crack  propagates  into  the  plate  normal  to  the  maximum  tangential  stress  near 
the  crack  tip.  If  This  maximum  stress  is  caused  to  shift  from  its 
original  direction,  the  direction  of  the  crack  propagation  would  be 
expected  to  change.  This  could  be  caused  by  the  stress  waves  generated 
by  the-  impact  and  is  probably  the  reason  for  A  not  following  a  straight 
line.  The  maximum  speed  of  the  crack  propagation  was  analytically 
determined  by  Yoffee  to  be  approximately  0.6  times  the  velocity  of  the 
shear  wave.  This  agrees  well  with  the  value  of  0.38  times  the  dilata- 
tional  wave  velocity  as  observed  by  KUppers  (1967)  in  experiments  with 
glass . 

It  is  apparent  that  the  stress  wave  associated  with  the  impact 
travels  through  the  plate  and  reflects  from  the  rear  surface,  inter¬ 
secting  the  crack  before  it  has  time  to  propagate  through  the  target. 
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It  is  computed  that  the  crack  reaches  the  first  point  of  branching  in 
approximately  3.5  n-sec  and  that  the  reflected  dilatational  wave  reaches 
the  sane  point  in  approdaately  3.4  j»-sec.  The  stresses  associated  with 
the  arrival  of  this  wave  will  casse  a  change  in  the  state  of  stress 
arotnd  the  propagating  crack  tip  and  is  probably  the  cause  of  the 
branching  which  creates  crack  C.  As  the  reflected  shear  wave  velocity 
is  less  than  that  of  the  dilatational  wave,  it  will  intersect  crack  A 
at  a  later  tine  and  would  be  expected  to  cuase  another  branching  such 
as  crack  B. 
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SLED  TEST  RESULTS 

A  rocket  sled  test  program  was  conducted  at  So llcatan  Air  Force  Base, 
Sew  Mexico,  to  detentice  tire  rain  erosion  resistance  of  slip-cast  fused 
sill -a  at  velocities  afccve  5000  ft/sec.  The  fabrication  techniques  used 
at  the  Georgia  Institute  of  Technology  to  slip-cast,  heat-treat,  and 
flaue-glaze  the  radoanes;  a  description  of  the  test  facilities;  and  the 
results  of  six  sled  tests  have  been  described  by  Walton  and  Karris  (1966). 

The  radome  from  run  Mo.  2,  Ho  Honan  Mo.  7RB1,  serial  Mo.  272-1  has 
been  examined,  and  some  observations  will  be  sade  concerning  the  inpact 
damage  that  it  sustained  in  the  40G-ft  rain  field-  The  rain  intensity 
was  reported  as  2.5  inches  per  hour  and  the  average  drop  size  was  said 
to  be  between  1.5- and  2.0-na  diameter.  The  average  velocity  was  about 
5300  ft/sec. 

Figure  6.1  shows  the  radome  on  the  sled  after  the  run.  A  view  of 
the  radome  with  its  metal  tip  broken  off  and  a  closer  view  of  some  of 
the  craters  are  shown  in  Figure  6.2.  Its  dimensions  as  well  as  the 
locations  of  the  centers  of  the  photographs  of  Figures  6.4,  6.5,  and 
6.6  are  given  in  Figure  6.3. 

The  photographs  of  Figure  6.4  were  taken  around  the  circumference 
at  45-degree  intervals  at  a  distance  of  11.6  inches  from  the  tip.  By 
counting  all  craters  that  could  be  detected, it  was  found  that  they 
averaged  27.1  per  square  inch.  Those  of  Figure  6.5  were  8.3  inches 
from  the  tip  and  averaged  40.4  craters  per  square  inch.  At  a  distance 
of  3.0  inches  from  the  tip, the  photographs  overlapped  and  they  are 
shown  in  Figure  6.6.  These  averaged  25.1  craters  per  square  inch.  At 
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tils  location  there  appeared  to  be  acre  erosion  from  the  inpacts  than 
craters.  Based  on  all  of  these  photographs,  it  was  estimated  that  far 
the  entire  radon e  there  sas  an  overall  average  of  31  craters  per  square 
inch.  There  is  probably  some  reason  as  to  why  the  crater  density  is 
greater  near  the  Middle  of  the  radoae  than  near  either  the  tip  or  the 
base*  but  we  have  no  explanation  at  this  tine.  The  radone  had  been 
f lane-glazed  by  the  use  of  a  plasna  jet  to  fuse  the  sjrface  and  tc 
provide  a  coating  of  non-porous  fused  silica  that  served  to  seal  the 
surface.  Walton  and  Harris  reported  the  thickness  of  the  fused  layer 
varied  fro*  0.03  to  0.04  inch  near  the  tip  to  less  than  0.005  inch  near 
the  base.  They  also  reported  that  although  each  drop  appeared  to  cause 
some  surface  damage  it  did  not  appear  to  extend  significantly  beneath 
the  glazed  surface.  This  is  in  agreement  with  our  findings  as  we 
measured  only  a  few  depths  greater  than  0.005  inch. 

Because  the  impact  velocities  of  all  drops,  both  large  and  small, 
are  essentially  the  same,  various  authorities  argue  that  the  impact 
pressures  as  given  by  the  "water-hammer’'  equation,  P  =  pCV,  will  be 
the  same  for  all  drops  where  p  is  the  density  of  the  liquid,  C  the 
velocity  of  the  compression  wave  in  the  liquid,  and  V  the  component 
of  the  impact  velocity  normal  to  the  surface.  Brunton  (1967)  explains 
that  one  of  the  reasons  why  a  large  drop  is  more  damaging  than  a  small 
one  at  the  same  velocity  is  that  with  the  large  drop  the  maximum  pres¬ 
sure  acts  for  a  longer  period  of  time.  This  would  also  result  in  longer 
stress  waves  being  propagated  into  the  target,  affecting  the  internal 
cracking  and  spall  damage.  The  areas  of  the  Craters  resulting  from  the 
larger  drops  would  obviously  be  greater  than  those  caused  by  small  drops. 
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Tbene  still  rewainr  -die  question  as  to  whether  or  not  the  "water  h-ti" 
pressure  is  die  nasi—. 

Crater  Characteristics 

the  depths  of  some  of  the  craters  were  neasured  and  depth  contours 
were  plotted.  The  purpose  of  this  is  to  learn  more  of  the  nature  of  the 
fractures  and  to  determine  the  extent  to  which  drop  impacts  can  be  simu¬ 
lated  by  the  use  of  the  water  jet  accelerator  and  by  solid  projectile 
iapacts.  Cross-sections  of  a  typical  crater  near  the  base  of  the  radoae 
are  shown  in  Figure  6.7.  The  depths  were  neasured  from  the  curved  sur¬ 
face  of  the  radone.  The  first  profiles  were  taken  outside  the  crater 
and  indicate  die  roughness  or  irregularities  of  the  surface.  As  it  is 
impossible  to  recognize  the  actual  crater  bomdaries  from  the  measure¬ 
ments  (although  they  can  be  seen  visually)  because  of  the  surface  rough¬ 
ness,  the  depths  shown  in  Figure  6.8  start  with  the  0.001-inch  contour. 
Comparing  the  profiles  and  depth  contours  of  this  crater  with  those 
produced  by  the  jet-produced  cavity  shown  in  Figures  5.9  and  5.10  several 
similarities  and  differences  will  be  noted.  In  both  cases  there  is  a 
central  plateau  surrounded  by  areas  of  greater  depth  containing  several 
rather  deep  pits.  The  primary  differences  are  that  the  radoroe  crater 
is  more  irregular  in  shape  (although  its  actual  boundary  is  not  shown) 
and  the  crater  surface  is  rougher.  The  jet-produced  crater  is  somewhat 
larger  and  deeper  than  that  formed  by  the  drop.  Also,  the  central 
plateau  of  the  first  was  the  original  surface  which  had  been  painted 
black  and  it  is  apparent  that  at  least  some  of  the  original  radome 
surface  had  been  removed  from  this  region  in  the  second.  The  differences 
are  simple  to  explain;  in  fact,  the  similarities  are  more  than  would  have 


been  expected. 
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Althoqgb  both  target*  were  of  slip-cart  fared  silica,  the  surfaces 
were  entirely  different.  As  stated,  the  redone  had  been  provided  with 
a  coating  of  nae-porons  fused  silica  which  varied  from  0.005  to  0.030 
inch  thick,  talton  (1966)  observed  that  race  the  inpact  pressure  was 
sufficient  to  cause  danny  in  the  glased  surface  that  there  was  a  ten¬ 
dency  to  chip  out  a  portion  of  the  glaze  extending  beyond  the  actual 
area  of  inpact.  In  an  unglazed  surface,  be  reported,  the  dsnage  was 
restricted  to  the  area  of  inpact,  iar  a  conparisoa  of  the  dsnage  re¬ 
sulting  from  jet  inpacts  and  drop  inpacts,  it  is  necessary  to  have  two 
identical  Materials.  Fran  these  results,  however,  it  appears  that  the 
drop  inpact  can  be  s  isolated. 

Another  explanation  of  the  difference  between  the  two  craters  is 
that  the  jet  inpact  was  nomal  to  the  target's  surface  at  a  velocity 
of  about  2000  ft/sec  and  chat  of  the  drop  was  at  an  angle  cf  about  15° 
at  a  speed  of  approximately  5000  ft/sec,  giving  a  normal  velocity  com¬ 
ponent  of  approximately  1300  ft/sec.  This  difference  may  explain  the 
difference  between  naxinun  crater  depth  of  0.008  inch  in  the  jet- 
created  crater  and  only  0.005  inch  in  the  drop-generated . 

The  direction  of  impact  seems  to  sake  little  difference  in  the 
crater's  general  appearance.  It  appears  that  the  simulation  of  a  high- 
velocity  impact  at  a  small  angle  with  the  surface  can  be  accomplished 
with  a  lower  velocity  impact  at  a  different  angle  as  long  as  the  normal 
velocity  components  are  the  same.  For  example,  the  velocity  of  5000  ft/ 
sec  at  15°  with  the  surface  could  probably  be  simulated  with  a  velocity 
of  2000  ft/sec  at  an  angle  of  40°. 

The  profiles  and  crater  contours  of  another  crater  located  near 
the  midpoint  of  the  radome  surface  are  given  in  Figures  6.9  and  6.10. 
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Although  this  crater  area  is  less  than  the  previous  one,  its  depth  is 
somewhat  greater. 


The  Rain  Field 

It  has  been  suggested  that  soae  of  the  drops  say  be  deflected  away 
froei  the  radome's  surface  by  the  hypersonic  shock  layer.  The  simulated 
rain  rate  for  the  Hollosan  tests  was  adjusted  to  2.5  inches  per  hour. 

The  drop  size  distribution  for  this  facility  as  given  by  Martens en  (1970) 
agrees  reasonably  well  with  the  one  for  natural  rain  given  by  Fyall  (1970). 
Based  upon  this  distribution  and  a  terminal  drop  velocity  of  6  m/sec,  also 
from  Fyall  (1970),  which  gives  a  concentration  of  2.94  gm/m3,  the  follow¬ 
ing  distribution  of  drop  sizes  was  computed: 


i  Diameter  (mm) 

Number  in  In 

0.25  - 

0.75 

672 

0.75  - 

1.25 

305 

1.25  - 

1.75 

241 

1.75  - 

2.25 

145 

2.25  - 

2.75 

74 

2.75  - 

3.25 

32 

3.25  - 

3.75 

15 

3.75  - 

4.25 

5 

4.25  - 

4.75 

2 

4.75  - 

5.25 

1 

1492  drops /ra3 

The  base  area  of  the  radome  times  the  length  of  the  rainfield  gives  a 

o 

volume  of  3.54  m  or  5147  drops  in  its  path.  Considering  the  radome  as 
conical,  its  surface  area  is  155  in.  .  This  gives  an  average  of  33  impacts 
per  square  inch  if  all  drops  in  the  radome  path  were  intercepted.  As  the 
number  of  craters  was  estimated  to  be  31  per  square  inch,  it  is  concluded 
that  practically  all  drops  in  the  radome' s  path  produced  craters.  This 
was  also  the  conclusion  reached  by  Walton  (1967),  who  stated  that  an 
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exantination  of  the  radome  revealed  that  the  number  of  drops  impacted 
corresponded  to  the  number  of  drops  larger  than  0.5  m.  It  is  also  in 
agreement  with  the  analytical  conclusions  reached  in  Section  II  of  this 
report  which  indicated  that  for  all  velocities  of  this  radome  between 
2000  and  10,000  ft/sec,  the  maximum  deflection  of  the  raindrops  would 
be  less  than  one-half  degree. 
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SECTION  VII. 

PLANNED  FUTURE  INVESTIGATIONS 

The  research  described  in  this  report  will  be  extended  in  several 
ways.  These  are  described  below. 

The  Motion  of  Raindrops  in  a  Hypersonic  Shock  Layer 

» 

First,  calculations  will  be  carried  out  for  blunt-nosed  (spherical) 
vehicle  shapes.  Since  the  shock  wave  is  detached  from  the  body  surface 
in  such  cases,  one  would  expect  that  the  effectiveness  of  the  inter- 
phase  drag  force  in  slowing  down  the  drops  before  they  strike  the 
vehicle  surface  would  be  greater  for  a  spherical  shock  layer  than  for 
a  conical  shock  layer.  This,  in  turn,  would  lead  to  less  likelihood 
of  damage  to  the  radome.  The  computations  described  above  will  deter¬ 
mine  what,  if  any,  lessening  of  susceptibility  to  rain  erosion  damage 
would  be  produced  by  nose  blunting  so  that  this  could  be  weighed 
against  other  factors  in  missile  design. 

Second,  the  interphase  drag  and  heat  transfer  laws  will  be  modified 
to  apply  to  the  case  of  moderate  void  fractions  (the  void  fraction  is 
the  ratio  of  the  air  volume  to  the  total  mixture  volume).  All  previous 
work  has  been  restricted  to  the  case  of  large  void  fractions.  The  sen¬ 
sitivity  of  the  solution  variables  to  the  value  of  the  void  fraction 
will  be  found  and,  if  necessary,  the  appropriate  correction  factors  to 
the  large- void- fraction  solutions  will  be  determined. 

Third,  extensive  calculations  will  be  carried  out  to  determine  the 
sensitivity  of  the  results  to  the  distribution  of  drop  sizes  present 
within  the  mixture.  Previous  work  has  assumed  that  the  exact  size  dis¬ 
tribution  is  unimportant  and  what  the  real  rain  field  can  be  replaced, 
for  the  purposes  of  analysis,  with  a  fictitious  rain  field  in  which  all 
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drops  have  the  sane  size.  This  effective  size  is  usually  taken  to  be 
the  average  size  of  the  actual  drops.  The  limits  of  validity  of  this 
assumption  will  be  found  and,  if  necessary,  correction  factors  will  be 
determined. 

Fourth,  an  attempt  will  be  Bade  to  obtain  approximate  closed-form 
solutions  for  conical  and  spherical  shock  layers.  Such  formulas  would 
allow  rapid  estimation  of  the  solution  variables  for  design  purposes. 

The  trends  of  the  numerical  data  obtained  thus  far  suggest  several 
possible  approximations  that  would  simplify  the  basic  differential 
equations  enough  to  make  possible  analytical  solutions.  The  limits 
of  validity  of  the  approximate  solutions  will  be  determined  by  compari¬ 
son  with  the  exact  numerical  results. 

Cratering  and  Erosion  Mechanism 

Research  performed  under  the  present  contract  has  shown  that  the 
jet-impingement  model  is  a  very  useful  tool  in  the  investigation  of 
high-velocity  impacts  of  drops  against  solid  targets.  This  model  will 
be  improved  as  described  below. 

An  important  ingredient  of  the  model  mentioned  above  is  the  distri¬ 
bution  of  the  pressure  applied  to  the  target  by  the  impinging  drop. 

Past  work  has  employed  pressure  distributions  obtained  by  fitting  curves 
to  experimented  data  found  from  observations  of  jets  impinging  on  rigid 
surfaces.  A  better  pressure  distribution  will  be  found  by  extending 
the  work  of  Huang  (1971)  who  investigated  the  two-dimensional  low-speed 
impact  of  drops  of  various  shapes  on  a  rigid  plane  by  numerical  solution 
of  the  hydrodynamic  partial  differential  equations  governing  the  motion 
of  fluid  i:i  the  drop.  Huang's  method  will  be  applied  to  both  normal  and 
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oblique  high-speed  impacts  on  elastic  surfaces.  The  thermal  energy 
equation  will  be  solved  together  with  the  hydrodynamic  equations. 

This  will  allow  the  calculation  of  heat  transferred  to  the  target. 

This  is  expected  to  be  important  for  high-speed  impacts.  The  boundary 
conditions  at  the  target  surface  will  be  modified  by  a  procedure  simi¬ 
lar  to  that  used  by  Lee  and  Cheng  (1971)  to  account  for  the  elasticity 
of  the  target.  The  pressure  distributions  thus  obtained -will  then  be 
used  in  the  jet- impingement  model  to  obtain  further  information  about 
crater  depths  and  sizes  and  the  influence  of  various  factors  on  these 
quantities . 

Fractures  Resulting  from  Stress  Waves  Generated  by  Liquid  Impact 

Equations  have  been  derived  and  computer  programs  have  been  written 
that  will  enable  values  of  stress  in  an  impacted  target  to  be  computed 
if  the  pressure-time  relationships  on  the  surface  are  known.  The  pres¬ 
sures  created  by  both  the  drop  impact  and  the  lateral  jetting  action 
ore  considered  in  the  computation  of  the  stresses  resulting  from  the 
generated  stress  waves.  The  reflections  of  these  waves  from  the  inside 
surface  of  the  radome  are  also  taken  into  account.  This  program  will 
also  compute  the  stresses  resulting  from  any  number  of  multiple  impacts 
at  any  time  intervals  and  on  both  normal  and  inclined  surfaces. 

When  a  rain  field  is  described  (i.e.,  inches  per  hour,  duration, 
and  drop  size  distribution)  and  the  radome  profile,  thickness,  and 
velocity  are  specified,  it  will  be  possible  to  completely  determine 
the  probable  damage.  Before  this  can  be  done,  however,  the  items  men¬ 
tioned  in  the  previous  paragraphs  such  as  the  size  and  motion  of  drops 
behind  the  shock  wave,  the  effects  of  aerodynamic  heating,  the  distri¬ 
bution  and  duration  of  pressure  applied  to  the  target,  and  the  fracture 


It  appears  that  perhaps  the  greatest  insight  into  the  mechanics  of 

« 

rain  erosion  may  be;  gained  from  experimental  and  semi- empirical  investi¬ 
gations.  Experimental  results  are  certainly  necessary  to  establish  the 
validity  or  limitations  of 1  theoretical  analyses.  Although  only  prelim¬ 
inary  experimental  work  has  been  done  tp  date,  the  results  are  in  qual- 

•  }  '  * 

itative  agreement  with  the  theory  of  fracture  that  has  been  proposed. 

•  i 

This  theory  does  not  agree  with  those,  proposed  by  other  investigators , 

!  I 

but  it  is  the  only  one  that  ,explains  all  types  of  fractures  that  hlave 
been  observed  to  result  from  liquid  impact.  Further  Work  will  be  per¬ 
formed  to  give  quantitative  results.  Although  emphasis  is  to  be  placed 

1 

upon  damage  to, ceramic  type  radoraa  materials,  Lucite  is  being  used  as 
the  standard  material  in  some  Of  these  studies  because  its  properties 

•  i 

are  well  known,  because  of  its  transparency  and  easy  determination  of 

internal  damage,  !and  because  its  surfaces  can  be  artificially  roughened 

! 

to  determine  the  effect  of  surface  conditions  on  the  lateral  jetting 

I 

I  *  •  i 

and  resulting  surface  pressures.  Once  the  distribution  and  duration  of 
the  pressures  on  the  target  are  determined  for  various  parameters,  these 

can  be  applied  to  other  materials  to  determine  the  damages  that  would  be 

•  } 

caused. 

A  liquid  jet  accelerator  has1  been  constructed  to  simulate  raindrop 
impact.  The  first  shots  gave  velocities  of  Approximately  2000  ft/sec. 

It  is  believed  that  by  making  modifications  in  design;  the  velocities  can 
be  at  least  doubled.  It  will  also  be  possible  to  accelerate  two  drops, 
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simltaneously  so  that  the  effects  of  Multiple  iaoact  can  be  determined 
and  compared  with  the  theoretical  results.  By  the  use  cf  other  liquids 
such  as  mercury  it  will  be  possible  to  simulate  higher  impact  velocities. 

A  few  liquid  impacts  of  slip-cast  fused  silica  targets  have  beat 
made.  The  craters  seen  to  be  very  similar  to  those  formed  during  rocket 
sled  tests. 

A  profilometer  has  been  constructed  with  which  measurements  of  the 
craters  can  be  made.  A  series  of  experiments  will  be  performed  from 
which  data  the  drop  velocities  and  sizes  can  be  correlated  with  crater 
characteristics  such  as  volume  and  depth  for  various  angles  of  impact. 
These  experiments  will  also  include  the  effects  of  varying  the  temper¬ 
ature  of  the  targets.  The  effects  of  surface  roughness  will  also  be 
investigated. 

It  seems  that  the  effects  of  the  reflected  stress  waves  from  the 
inner  surface  is  very  important.  Spallation  may  present  a  serious 
problem  if  thin-walled  radomes  are  used.  Detachment  of  radome  material 
from  the  inner  surface  at  high  velocities  would  destroy  equipment. 

Water  impact  experiments  shew  that  the  rear  surface  of  the  target  often 
suffers  a  greater  damage  than  the  front  surface.  Experiments  have  also 
demonstrated  the  importance  of  having  a  smooth  and  uniform  rear  surface 
as  small  defects  such  as  scratches  and  voids  greatly  weaken  the  target. 

Another  phase  of  the  proposed  experimental  program  is  to  impact  the 
targets  with  small  solid  particles.  There  are  it  least  three  reasons  for 
this:  (1)  In  flight,  aircraft  and  missiles  encounter  foreign  bodies 
other  than  water  drops;  (2)  To  aid  in  a  more  complete  understanding  of 
the  mechanics  of  the  collision  process;  (3)  To  examine  the  merits  and 
limitations  of  simulation  techniques  in  which  solid  particles  are 
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substituted  for  water  drops.  The  last  of  these  is  perhaps  the  principal 
reason  for  including  these  experiments  in  this  program.  A  light-gas  gun 
with  the  capability  of  accelerating  drop-size  plastic  projectiles  to 
velocities  as  high  as  40*000  ft/sec  has  been  constructed  at  Tennessee 
Tech.  Lead  and  polyethylene  spheres  have  been  used  to  simulate  rain 
erosion.  Plasticine  (modeling  clay)  is  said  to  most  closely  resemble 
water  flow  characteristics  on  impact. 

A  ballistic  pendulum  will  be  used  to  determine  the  energy  and 
momentum  transfer  to  the  target.  The  momentum  multiplication  factor  is 
a  very  important  value  to  be  determined  for  both  normal  and  oblique 
impacts.  Other  types  of  instrumentation  will  include  the  use  of  strain 
and  thermal  gages,  interferometry,  photoelasticity,  etc.  as  needed  to 
identify  the  important  parameters  that  are  critical  to  the  radome's 
performance  in  a  rain  environment. 
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FIGURE  2.21.  SURFACE  DISTRIBUTIONS  OF  DROP-PHASE  NORMAL  VELOCITY 
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FIGURE  2.23.  SURFACE  DISTRIBUTIONS  OF  DROP-PHASE  NORMAL  VELOCITY 
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FIGURE  3.5.  SHEAR  FACTOR  EFFECTS  IN  STRONG  PENETRATION  (V*80Q0  CM/SEC) 
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FIGURE  3.8.  IMPACT  OF  WATER  JET  ON  CARBON  TETRACHLORIDE  SURFACE 
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FIGURE  3.11.  IMPACT  OF  MERCURY  DROPS  ON  2024-0  ALUMINUM 


MAXIMUM  CAVITY  DEPTH 


MAXIMUM  CAVITY  DEPTH 


0 


.2# 


.24 


.22 


.20 


to 


.  14 

.14 


.12 


.10 


.08 


.08 


.04 


.02 


Colculoted  Results 

Experimental  Data  off 
Engel  (I960) 


1  .  ,  I _ l _ i _ I _ i _ t _ i_- _ i _ i 

2  3  4  5  6  7  8  9  10  11*10* 


IMPINGEMENT  VELOCITY  (  CM/SEC 


FIGURE  3.13.  IMPACT  OF  MERCURY  DROPS  ON  COPPER 
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FIGURE  3.15.  IMPACT  OF  MERCURY  DROPS  AGAINST  STEEL 


FIGURE  4.1  FRACTURES  CAUSED  BY  REFLECTED  STRESS  WAVES 


FIGURE  4.2  PRINCIPAL  STRESSES  PROOUCED  BY  IMPACT 
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FIGURE  5.1  WATER  JET  ACCELERATOR 


FIGURE  5.2  IMPACT  RANGE 


FIGURE  5.2  -  CONTINUED 


FIGURE  b.3  HIGH-SPEED  WATER  JET 


FIGURE  5.4  WATER  IMPACT  DAMAGE  TO  LUCITE  TARGETS 
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FIGURE  5.7  PROFILE  AMD  DEPTH  CONTOURS  OF  CRATER  IN  LUCITE 
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FIGURE  5.8  WATER  IMPACT  DAMAGE  TO  FUSED  SILICA  TARGETS 
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FIGURE  6.3  DIMENSIONS  OF  RADOME 


FIGURE  6.4  PHOTOGRAPHS  OF  CRATERS  NEAR  BASE 
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FIGURE  6.7  CRATER  CROSS  -  SECTIONS 
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APPEXDIX  1 


K.  V.  Bennett,  D.  E.  Re:  ioids,  and  J.  C.  Saelcer 


The  purpose  of  this  appendix  is  to  present  a  closed-form  approximate 
solution  for  the  drop-phase  flow  field  on  thin  wedges  and  cones.  Such 
solutions  are  useful  in  that  they  allow  designers  to  make  rapid  estimates 
of  the  flow*  variables  without  having  to  engage  in  extensive  nuierical 
calculations .  The  notation  is  the  sane  as  in  section  II  of  the  body  of 
this  report  unless  otherwise  stated.  For  simplicity  in  the  analytical 
work  attention  is  restricted  to  -the  linear  drag  law  and  tangential  and 
normal  coordinates,  x  and  y  respectively,  are  used. 

In  dimensionless  form  the  governing  equations  are 


«V  +  n)Vp’’s  +  (<V  *  *  0 

F  F  ,  +  G  F  ,  =  o,(l  -  F_) 

p  p’s  p  p  n  lv  P 


Vp's  *  GpVn  *  -»!<«”/=>  *  V 
Vp's  *  GpHp-„  *  “l(Hs  *  Hp> 

where  (Ala)  results  from  balance  of  mass,  (Alb)  from  balance  of  tangential 
linear  momentum,  (Ale)  from  balance  of  normal  linear  momentum,  (Aid)  from 
balance  of  energy,  and  the  dimensionless  variables  are  defined  as  in  (2.9) 
and  (2.10).  Defining  a  stream  function  such  that 

*P's  *  -<ebstn)Vp>  ♦p'n  *  <Vtn)Vp 

and  converting  the  governing  equations  to  independent  variables  ip  and 

z  =  n  one  obtains 

G  F  ,  =  a, (1  -  F  ) 

p  p’z  1  p 

VVz  =  -ai«3z/s)  +  Gp) 

^p^p’z  =  a2^Hs  “ 


A-2 


Gs  j  —  F 
P  z  P 


«p  =  ♦  z,\s’t  > 

p 

Using  (A2)  evaluated  cn  the  shock  a  relation  between  the  origin  of  each 


(A3) 


streamline  z  and  the  stream  function  can  be  found.  This  was  used  to 
s 

rewrite  (A3e)  as 


<L  =  -Os/(es-eb))1*i(zs/(ebs+z))i(i/G  s,z  )  (A4) 

‘  ^  s 

The  appropriate  boundary  conditions  on  the  shock  z  =  zs  are 

V1’  VV  S'1'  S'1'  s  =  v<vv  (A5) 

It  was  first  attempted  to  solve  (A3)  and  (A4)  by  a  perturbation 
method  for  small  a^.  The  first  two  approximations  we~>e  found  but  the 
labor  .involved  in  computing  further  approximations  was  so  extensive  that 
this  was  not  done.  The  second  approximation  was  found  to  be  valid  only 
for  «  1.  Thus  the  perturbation  scheme  was  abandoned.  Instead  it  was 
decided  to  linearize  the  nonlinear  terms  in  (A3)  using  the  values  of  Gp 
and  s  on  the  shock.  When  this  is  done  the  linearized  equations  have  the 
form 


-etGp’z  =  -a1(j(es-eb)(z/zs)  +  Gp) 

-ebHp’z  =  “2(Hs  -  V 
-0bs’z  =  Fp 

The  solutions  to  (A6)  satisfying  the  boundary  conditions  (A5)  are  easily 
found  to  be 


( A6) 


Gp  =  -0bexp(-(a1/eb)(zs~z))  +  j((es~eb)/zs)((zs 
+  (9b/a1))exp(-(a1/6b)(zs-z))  -  (z  +  (6^^))) 
Hp  =  Hs  -  (Hs-l)exp(-(a2/eb)U-0) 
s  =  (zs/(es-eb))  +  ((zs-z)/eb) 


( A7 ) 
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Using  (A7)  and 

*  W.V 

the  dimensionless  density  can  be  easily  evaluated  from  (A*»). 

Figures  A.l  and  A. 2  shew  comparisons  of  the  approximate  results 
with  the  numerical  solutions  obtained  for  slender  wedges  and  cones  in 
section  II.  The  density  ratio  is  e  =  0.1  for  this  data.  The  quantities 
chosen  for  comparison  are  the  normal  velocity  and  the  density  of  the 
drop  phase  evaluated  at  the  vehicle  surface.  It  can  be  seen  that  for 
small  and  moderate  values  of  cs^  the  approximate  solution  in  good 
agreement  with  the  numerical  results  (for  the  wedge  the  numerical  results 
are  identical  to  the  exact  analytical  solution  given  in  section  II).  For 
=  10  the  approximate  solution  begins  to  overestimate  the  impact  veloc¬ 
ity  and  underestimate  the  density.  This  could  probably  be  corrected  by 
modifying  the  linearizing  factors  but  this  was  not  attempted  in  the 
present  work. 

The  agreement  between  the  exact  and  approximate  solutions  appears 
to  be  satisfactory  for  small  and  moderate  values  of  the  drag  parameter 
a^.  It  appears  that  a  similar  procedure  could  also  be  used  to  obtain 
approximate  results  for  general  shock  layer  problems. 
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